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ABSTRACT 

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint 
Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars 
within 67 pc (10', or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes 
consistent with membership, down to a magnitude limit of r = 21.7. Based on photometric, kinematic, 
and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far 
out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, 
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we determine a velocity dispersion of ?>.T_ii km s . The mass within the half-light radius is 5. 
10^ M0. The stellar kinematics of Segue 1 require very high mass-to- light ratios unless the system is 
far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed 
toward implausibly short periods. With a total luminosity less than that of a single bright red giant 
and a V-band mass-to-hght ratio of 3400 M0/L0, Segue 1 is the darkest galaxy currently known. We 
critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic 
samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < —3) 
of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively 
that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that 
contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest estimated 
dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter 
physics and galaxy formation on small scales. 

Subject headings: dark matter — galaxies: dwarf — galaxies: kinematics and dynamics — galaxies: 
individual (Segue 1) — Local Group 



INTRODUCTION 

Survey (SDSS) has 



The Sloan Digital Sky 
tremendously successful in revealing new 
Way dwarf g alaxie s over the past five years (e.g., 
Willman et al.il200a IZucker et al.ll2006t IBelokurov et alj 



been 
Milky 



2007at IWalsh et all 120071: IBelokurov et all 120101) . How- 
ever, its limited depth and sky coverage, along with the 
difficulty of obtaining spectroscopic followup observa- 
tions, still leave us with an incomplete understanding 
of the Milky Way's satellite population. In particular. 
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key parameters such as the luminosity function, mass 
function, radial distribution, and total number of 
satellites depend extremely sensitively on the properties 
of th e few least luminous dwarfs (e.g., iTollerud et al.l 
l2008f ). which are not yet well-determined. Since the 
least luminous dwarfs are the closest and densest 
known dark matter halos to the Milky Way, these 
same objects represent critical target s for indirec t dark 



matt er detectio n experiments (e.g., 'Ba ltz et al.l 
Evans. Ferrer 



200C 



2007 



& Sarkar 2004; Colafrancc sco et al. 
Strigari et al. 2008b; Kuhlen, Dicmand, & Madau 200? 
mann ct al. 2009 : Fieri et al. 2009: Martinez et al.l 



Brin 



20091 ) and for placing limits on the phase space density 



of d ar k matte r particles (e.g. 
2000t iDalcanto n & H oganI [200lt 



Simon fc Gehal l2007l : iStrigari et al.l l2008bl : iGeha et al.l 
200il ). However, as the closest known satellites to the 
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Milky Way, they are also the most susceptible to tidal 
forces and other observational systematics. 

Because of the extreme lack of bright stars in these 
systems, most of the fainte st dwarfs su ch as Will- 
man 1 (iWillman et al.l I2005D. Bootes I I (jWalsh et al.l 
"2007), Segue 1 (Bclok urov et all l2007aD . and Segue 2 
(Bclokurov et al. 2009 ) remain relatively poorly char- 
acterized by observations; for example, the dynamical 
state of Willman 1 has still has not b een established 
(jMartin et al.ll2007Hw"illman et al.ll2010D . and the veloc- 
ity dispersion of Boo II is uncertain at the facto r of ^ 5 
level ( Koch et al. 2009). Similarly, although Geh a et all 
(|2009l hereafter G09) demonstrated that the kinemat- 
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ics of stars in Segue 1 clearly indicate that it is a dark 
matter-dominated object, other observations have sug- 
gested the possibility of tidal debris in the vicinity of 
Segue 1, as well a s potential contamination fro m the 
Sagittarius stream (jNieder ste-Ostholt et alJ[2"00 9t). 

More generally, the issues o f tidal disr uption (e.g., 
iPeharrubia. Navarro, fc McConi iachiel I2008D and binary 
stars (|McConnachie &: C6t6i2010l ) are the last remain- 
ing major questions to be settled regarding the nature 
of the faintest dwarfs. These objects promise cl ues to 
the extreme limits of gala xy formation (Gilmor e et al.l 
I2007t IStrigari et al.l [2008al ) and perhaps to the forma- 
tion of the first gala xies in the early universe (e.g., 
iBovill &: Ricottil [2009D . as well as offering insights into 
dark matter physics. However, these applications hinge 
on the assumption that the mass distribution of each sys- 
tem is accurately known. Current mass estimates assume 
dynamical equilibrium and that the observed kinematics 
are not being affected by Galactic tides or binary stars, 
but tests of those assumptions are obviously required in 
order to confirm that the dwarfs are bound, equilibrium 
systems. If instead the observed velocity dispersions of 
Segue 1, Willman 1, and others are being inflated either 
by the tidal influence of the Milky Way or the presence of 
binary stars in the kinematic samples, then they are un- 
likely to be useful probes of the behavior of dark matter 
on small scales. 

Correcting velocity dispersions for binaries, which are 
inevitably prese nt in any stellar s ystem, is relatively 
straightforward ([Minor et al.l l201Clf ). The only obser- 
vational requirement is that a significant subset of the 
sample have at least two velocity measurements with a 
separation of order 1 yr. Tidal effects, unfortunately, 
are more difficult to nail down. The only unambigu- 
ous signature of tidal interactions is the presence of 
tidal tails (e.g.. lToomre fc Toomrel[l972| ). Detecting such 
features in the ultra-faint dwarfs is extremely challeng- 
ing: the galaxies themselves ha ve central surface bright- 
nesses of 26 — 28 mag arcsec"^ ([Martin, de Jong, fc Rixl 
[2008), so any tidal debris would be at least several 
magnitudes fainter and likely below the SDSS detec- 
tion limit of ~ 30 mag arcsec"^. Deeper, wide- 
field photometric surveys of the ultra-faint dwarfs can 
reach surfac e brightnes ses as low as 32.5 mag arcsec"'^ 
(ISand et al.l[20b 9. 2010: IMufioz. Geha! fc Willmanl[20Tol: 
Ide Jong et al.l[201Q) . but such observations are not yet 
available for most of the dwarfs. 

In principle, spectroscopic studies can pinpoint the 
stars associated with an object and probe debris at lower 
surface densities than is possible photometrically. Spec- 
troscopic surveys also provide the only means of identify- 
ing tidal d ebris that is orient ed along the line of sight to 
an object ([Lokas et al.]l2008t iKlimentowski et al.l 120091 ). 
However, the currently available spectroscopic samples 
of less than 25 stars in the faintest dwarfs are not suffi- 
cient to determine to what extent tides may be affecting 
the kinematics. Much larger spectroscopic data sets are 
required to test for tidal effects. 

In this paper, we present a nearly complete spectro- 
scopic survey of Segue 1 that is aimed at obtaining re- 
peated velocity measurements of known members and 
searching for stars that have been tidally stripped from 
the system. We describe our modeling of the binary star 
population and the mass distribution in more detail in 



a companion paper ([Martinez et al.l [20101 hereafter Pa- 
per II), and a separate study examines the implications 
of our new m ass measurements for indirect detection of 
dark matter ([Essig et al.l I2010D . In Section [2| we de- 
scribe the survey and the data reduction. We identify 
Segue 1 member stars in Section [31 and then analyze 
their metallicities and velocities in Section [H In Sec- 
tion [5] we present our derivation of the intrinsic velocity 
dispersion of Segue 1 after correcting for the presence 
of binary stars in the sample (see Paper II for more de- 
tails), and in Section [51 we describe our detection of an 
unrelated tidal stream in the same part of the sky. We 
consider the implications of this data set for proposals 
that the kinematics of Segue 1 are affected by contami- 
nation and tidal disruption in Section [71 We discuss the 
utility of Segue 1 for placing constraints on the properties 
of dark matter in Section [51 In Section [31 we summarize 
our findings and conclude. 

2. EXPERIMENTAL DESIGN, OBSERVATIONS, 
AND DATA REDUCTION 

2.1. A Survey for Tidal Debris 

As a complement to ongoing deep, wide-field 
photometric surve ys of the ultra- faint dwarfs (e.g., 
iMuhoz et al.ll20iol ). we embarked upon a spectroscopic 
search for evidence of tidal stripping or extratidal stars. 
The ideal target for such a search would be a galaxy that 
(1) is nearby, to maximize the tidal forces it is currently 
experiencingl3 (2) is moving at a high velocity relative to 
the Milky Way, to minimize the degree of contamination 
by foreground stars, and (3) has a small angular size, to 
minimize the area that the survey needs to cover. Out of 
all the known Milky Way dwarf galaxies, the clear choice 
according to these criteria is Segue 1. At a distance of 
23 kpc from the Sun (28 kpc from the Galactic Center), 
Segue 1 is the closest dwarf galaxy other than Sagittarius, 
which of course is the prototype for a dwarf undergoing 
tidal disruption. Its heliocentric velocity of 207 km s~^ 
(the largest of the Milky Way satellites within 200 kpc) 
and relatively small velocity dispersion give Segue 1 the 
lowest expected surface density of Milky Way foreground 
stars within 3(7 of its mean velocit y (according to the Be- 
sangon model; IRobin et al.ll2003l) . Finally, if Segue 1 is 
not surrounded by a massive dark matter halo — and it 
can only host visible tidal features if no extended halo is 
present — its instantaneous Jac obi (tidal) rad i us bas ed 
on the stellar mass estimated bv iMartin et al.l ([2008D is 
^ 30 pc, or 4.5', which is an observationally feasible area 
to search. This calculation conservatively assumes that 
Segue 1 has never been closer to the Milky Way than 
it is now; if its orbital pericenter is less than 28 kpc, 
its baryon-only tidal radius would be even smaller. The 
properties of Segue 1 are summarized in Table [1] 

2.2. Target Selection 

To select targets for the survey, we focused on the 
area within ^ 15' (100 pc) of the center of Segue 1 as 
determined bv IMartin et al.l ([20080 . Guided by the 24 
member stars identified by G09, we tweaked the color of 
the (appropriately shifted and reddened) M92 isochrone 

^ If the object is too close to the pericenter of its orbit, though, 
then the extent of its tails (if they exist) would be minimized. 
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TABLE 1 
Summary of Properties of Segue 1 



Row 


Quantity 


Value 


(1) 


RA (J2000) 


10:07:03. 2±1. 7' 


(2) 


Dec (J2000) 


+16:04:25±15" 


(3) 


Distance (kpc) 


23 ±2 


(4) 


Mv,o 


1 t: + 0.6 

J-3_o.8 


(5) 


Lv,o (L©) 


340 


(6) 


e 


0.48«:f3 


(7) 


fiv,o (mag arcsec-2) 


27.6tS:? 


(8) 


^off (pc) 


29t^ 


(9) 


Vhe! (km s--^) 


208.5 ±0.9 


(10) 


Vgsr (kms 1) 


113.5 ±0.9 


(11) 


a (km s"-*^) 


3.7111 


(12) 


Mass (Mq) 


5.8^11 X 10^ 


(13) 


M/Lv (M0/L0) 


3400 


(14) 


Mean [Fe/H] 


-2.5 



Note. — Rows (1 
taken from the SDSS 
of [M artin ct al. (2008) 
iBelokur ov et ak (2007a). 
are derived in this paper. 



)-(2) and (4)-(8) are 

photometric analysis 

and row (3) from 

Values in rows (9)- (14) 



from IClem et alj (|2008D so that it passed through the 
center of the member sequence at all magnitudes (this 
adjustment was only needed for the subgiant branch and 
main sequence, not the red giant branch). The Segue 1 
main sequence appears to be slightly redder than that 
of M92, with the offset increasing toward fainter mag- 
nitudes (see Figure [ij. The full G09 member sample is 
located within 0.25 mag of the adjusted fiducial track 
(0.2 mag for r < 21 and 0.1 mag for r < 20). 

Using po sitions and magnitudes extracted from DR5 of 
the SDSS (jAdelman-McCarthv et~aIll2007D . we selected 
stars within a narrow range of colors around the adopted 
fiducial sequence. Red giants (r < 20) were required to 
be within 0.1 mag of the sequence, and the selection win- 
dow was gradually widened toward fainter magnitudes, 
reaching 0.235 mag at r = 21.7. Horizontal branch can- 
didates were allowed to be 0.2 mag away from the fiducial 
track. A small number of stars located near a metal-poor 
asymptotic giant branch isochrone from IGirardi et al.l 
11004) were also selected. Within 10' (67 pc) of Segue 1, 
we identified 112 stars lying within the color- magnitude 
selection box (down to a magnitude limit of r = 21.7) 
that we consider to be our primary target sample. Stars 
up to a factor of two farther away from the fiducial se- 
quences or within the primary color-magnitude selection 
region but at larger distances from Segue 1 were targeted 
with reduced priorities. We also included as many of the 
known member stars as possible on multiple slitmasks to 
obtain repeated velocity measurements for constraining 
the binary population. 

2.3. Observations 

We observed 12 new slitmasks, including at least 
one slit placed on each of the 112 candidate member 
stars (plus repeat observations of 18 of the 24 members 
from GOD), with the DEIMOS spectrograph (jFaber et al.l 
120031) on the Keck II telescope. The observations took 
place on the nights of 2009 February 18, 26, and 27 and 
2010 February 12 and 13. 

The spectrograph setup and observing procedures were 




-0.5 



2.0 



Fig. 1. — Photometric selection criteria for candidate Segue 1 
members. The red giant branch/main sequence selection region 
(shaded red) is based on the M92 isochrone of Clem et aljl (2008), 
adjusted slightly at magnitudes fainter than r = 20.65 so as to en- 
close all of the spectroscopically confirmed members from G09 (red 
line). The blue and green shaded r egions rep resent the horizontal 
branch (from the M13 isochrone oflClem e t al.ll2008l) and asymp- 
totic giant branch (from a Girardi et al. 20(3 theoretical isochrone 
at [Fe/H] = —1.7) selection boxes, respectively. The filled points 
are the 24 member stars identified by G09. 

identical to those described by SG07: we used the 
1200 i?/mm grating with an OG550 filter to cover the 
wavelength range 6500 — 9000 A at a spectral resolution 
of i? = 6000 (slit width of 0.7"). An internal quartz lamp 
and Kr, Ar, Ne, and Xe arc lamps were employed for flat- 
fielding and wavelength calibration, respectively. Total 
integration times for the science masks ranged from 10 
minutes for a mask targeting very bright stars to ^ 2 
hours for most of the masks aimed at fainter stars. The 
masks are summarized in Tabled Conditions during the 
observing nights were generally good, with seeing rang- 
ing from 0.7 — 1.0" and thin cirrus at times. 

2.4. Data Reduction 

As with the observations, data reduction followed the 
outline given in SC07. We used a modified version of the 
data reduction pipeline developed for the DEEP2 galaxy 
redshift survey. The additional improvements we made 
to the code since SG07 were to re-determine the wave- 
length solution for slits that were initially not fit well and 
to identify and extract serendipitously observed sources 
more robustly. The pipeline determines a wavelength so- 
lution for each slit, flat-fields the data, models and sub- 
tracts the sky emission, removes cosmic rays, coadds the 
individual frames, and then extracts the spectra. 

In the spring of 2009, the DEIMOS CCD array was 
experiencing an intermittent problem wherein one of the 
eight chips (corresponding to the blue half of the spectral 
range at one end of the field of view) would fail to read 
out completely on some exposures. In a few cases, the 
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TABLE 
Keck/DEIMOS 


2 

Slitmasks 






Mask 


a (J2000) 


5 (J2000) 


PA 


icxp 


MJD of 


# of slits 


% useful 


name 


(h m s) 


(° ' ") 


(deg) 


(sec) 


observation 




spectra 


Seguel-l'^ 


10 07 06.01 


16 02 56.1 


65.0 


5400 


54416.551 


59 


88% 


Seguel-2 


10 07 08.85 


16 04 51.9 


180.0 


7200 


54881.311 


65 


83% 


Seguel-3 


10 07 00.82 


16 06 59.6 


-57.0 


5400 


54881.416 


61 


92% 


Seguel-C 


10 06 39.29 


16 06 02.0 


171.0 


600 


54881.284 


38 


87% 


segtidel 


10 06 57.65 


16 10 20.5 


178.0 


7500 


54889.267 


48 


88% 


segtide2 


10 06 36.07 


16 07 39.6 


171.0 


7500 


54889.343 


48 


90% 


segtideS 


10 07 34.37 


16 10 28.9 


144.0 


7500 


54889.444 


47 


94% 


segtide4 


10 06 19.07 


16 01 34.7 


179.0 


7350 


54890.243 


50 


84% 


segtideS 


10 07 40.78 


15 56 22.1 


1.0 


7800 


54890.342 


44 


89% 


segtide6 


10 06 59.81 


15 56 50.8 


-78.0 


7200 


54890.444 


54 


93% 


segtide? 


10 07 14.20 


15 55 43.9 


-123.0 


2700 


54889.542 


49 


76% 


segtideS 


10 06 58.45 


16 03 37.4 


-134.0 


1200 


55240.355 


40 


90% 


segtide9 


10 07 07.83 


15 56 24.7 


-165.0 


1800 


55240.402 


48 


85% 



'^ Seguel-1 is the slit mask observed by G09. For consistency, we use the velocities reported 
in that paper rather than re-reducing the mask and measuring the velocities again. 



affected cfiip was completely blank, while in others only 
the top or bottom of the chip was absent and the remain- 
der of the pixels were saved normally. We dealt with this 
problem conservatively by excluding all of the data from 
the temperamental chip from our reduction any time it 
exhibited abnormal behavior. As a result, a fraction of 
our targets have lower S/N in the blue than would be 
expected from the total mask observing times, but since 
most of the key spectral features for our analysis are on 
the red side of the spectra, the overall impact is minor. 
After reducing the data, we used the custom IDL code 
described by SG07 to measure the radial velocity of each 
star. We first corrected each spectrum for velocity off- 
sets that could result from mis-centering of the star in 
the slit by cross-correlating the telluric absorption fea- 
tures against those of a t elluric standard s tar specially 
obtained for this purpose ()Sohn et al.ll2007l SG07). The 
spectra were then cross-correlated with a library of high 
S/N templates with well-known velocities obtained with 
DEIMOS in 2006 and 2007. The template library con- 
tains 15 stars ranging from spectral type F through M, 
mostly focusing on low-metallicity giants but also in- 
cluding representative examples of subgiants, horizontal 
branch stars, main sequence stars, and more metal-rich 
giants. We also fit 4 extragalactic templates, identifying 
a total of 69 background galaxies and quasars. The veloc- 
ity of each target spectrum is determined from the cross- 
correlation with the best-fitting template spectrum. We 
estimate velocity errors using the Monte Carlo technique 
presented in SG07: we add random noise to each spec- 
trum 1000 times and then redetermine its velocity in each 
iteration. We take the standard deviation of the Monte 
Carlo velocity distribution for each star as its measure- 
ment uncertainty. Previous analysis of a sample of stars 
observed multiple times with the same instrument config- 
uration and analysis software indicates that our velocity 
accuracy is limited by systematics at the 2.2 km s~^ level. 
We have now confirmed the magnitude of the systematic 
errors with a much larger sample of repeat measurements 
than were used by SG07. 



2.5. Repeat Measurements 



As discussed in Sections [T] and 12.21 one of the main 
goals of this study is to determine the effect of binary 
stars on the observed velocity dispersion of Segue 1. Ac- 
complishing this task requires making multiple velocity 
measurements of stars, verifying that the derived velocity 
uncertainties are accurate, and searching for individual 
binaries. The bulk of our analysis of the repeat obser- 
vations is presented in Paper II, but in Figure ^Bp, we il- 
lustrate the agreement between subsequent velocity mea- 
surements and the initial one for each star observed more 
than once. Figure [2j) shows the distribution of velocity 
differences for each pair of measurements, normalized by 
their uncertainties. The good match to a Gaussian dis- 
tribution for —2 < {v2 — 'V\)l \/ (y\ + erf < 2 indicates 
that the uncertainties are accurate, and the excess in the 
wings of the distribution provides evidence for velocity 
variability (see § l4.2l and Paper II). A total of 93 stars, in- 
cluding approximately half of the member sample, were 
observed at least twice during the course of our survey. 

2.6. Spectroscopic Completeness 

In ? I2.3[ we described obtaining a spectrum of each 
one of the 112 stars within our primary photometric se- 
lection region and not more than 67 pc from the center of 
Segue 1. We successfully measured velocities for 109 of 
these stars. One target star (SDSSJ100707.12+160022.4) 
did not have any identifiable features in its spec- 
trum, one spectrum (SDSSJ100700.75-hl60300.5) suf- 
fered from reduction difficulties (see Appendix), although 
it appears to be a member, and the remaining one 
(SDSSJ100733.12-H155736.7) was not detected in our 
data. The magnitude of this latter star should have 
made it easily visible in the exposures we obtained, but 
a check of the SDSS images showed no source at this po- 
sition, and the target was flagged in the SDSS catalog as 
a possible moving object. We conclude that this source 
was actually an asteroid that happened to have the right 
colors and position to be selected for our survey. After 
excluding it, our completeness is 98.2%. We also tar- 
geted a fraction of Segue 1 member candidates as far out 
as 16' (107 pc). Within 11', 12', and 13', our overaU 
completeness is 97.0%, 96.7%, and 92.3%, respectively. 
Only a few stars at larger radii were observed. 
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Fig. 2. — (a) Comparison of velocity measurements for the sample of stars observed at least twice. The inset shows a zoomed-in view 
of the velocity range near Segue 1. (fe) Histogram of velocity differences between repeat observations, normalized by their uncertainties 

((d2 — fi)/i/a-^ + ctI). The red curve is a Gaussian with unit dispersion, which the data should follow if the measurement uncertainties 

are correct. The larger than expected number of stars in the wings of the distribution is caused by the presence of binaries and RR Lyrae 
variables in the sample. 



3. DEFINING THE SEGUE 1 MEMBER SAMPLE 

In total (including the observations of G09) we ob- 
tained 522 good spectra of 393 individual stars, of which 
167 were classified as high-priority member candidates 
(109 within 10' of the center of Segue 1 and 58 at larger 
radii) according to the criteria described in ii l2.2l We 
present all of our velocity and Ca triplet equivalent width 
measurements in Table 3. Our repeat velocity measure- 
ments for 93 of these stars span a maximum time baseline 
of 2.25 years. We use a variety of techniques to identify 
Segue 1 member stars in this data set. The primary data 
available for distinguishing members from non-members 
are: color/magnitude, velocity, metallicity (either in the 
form of [Fe/H] estimated from the Ca triplet lines or sim- 
ply the raw Ca triplet equivalent width), spatial position, 
and the strength of the Na I A8190 A doublet^ 

3.1. Methods of Identifying Segue 1 Members 

We consider three different membership selection tech- 
niques: a subjective, star-by-star selection using velocity, 
metallicity, color, magnitude, and the spectrum itself, a 
slightly modified ver sion of the algorithm introduced by 
iWalker et al.l (|2009bl ). and a new Bayesian approach pre- 
sented in Paper II. For the remainder of the paper, we 
adopt the third method as defining our primary sample 
except where explicitly noted. 

^^ We note, however, that the Na lines only function as a way to 
distinguish dwarfs from giants for very red stars (V — I > 2), signif- 
icantly redder than anything we expe ct to find i n a low-luminosity, 
low-metallicity system like Segue 1 (iGilbert et al. 2006). Also, be- 
cause of the tiny number of giants in Segue 1, the majority of the 
member stars are actually on the main sequence rather than the 
giant branch. Therefore, the strength of the Na doublet primarily 
serves to eliminate stars that were already obvious non-members 
from consideration. 



We first select member stars using the parameters 
listed above. Candidate members must have colors and 
magnitudes consistent with the photometric selection re- 
gion displayed in Figure [T] and described in Ji l2.21 must 
have velocities within a generous ~ 4cr window around 
Segue 1 based on the systemic velocity and velocity dis- 
persion measured by G09, and should have low metallic- 
ities. We can then iteratively refine the member selec- 
tion by examining the stars near the boundaries of the 
selection region more carefully. With this process, we 
classify 65 stars as definite members, 6 additional stars 
as probable members, and 5 more as likely (but not cer- 
tain) non-members. The remaining stars are clearly not 
members. The sample selected in this way is displayed 
in Figure [3] 

While the flexibility afforded by this subjective ap- 
proach is useful and allows all available information 
to be taken into account, a rigorous and objective 
method is als o desirable to avoid the possibility of bias. 
IWalker et al.l (|2009b) recently developed such a statis- 
tical algorithm to separate two potentially overlapping 
populations. This technique, known as expectation max- 
imization (EM), allows one to estimate the parameters of 
a distribution in the presence of contamination, with spe- 
cific application to the case of identifying dwarf galaxy 
member stars against a Milky Way foreground. Given 
a set of velocities, radial distances, and a third pa- 
rameter in which dwarf galaxy stars and Milky Way 
stars are distributed differently, the EM algorithm re- 
lies on the distinct distributions of the two populations 
in each property to iteratively assign membership prob- 
abilities to each sta r until it converges on a solution. 
IWalker et al.l ()2009bD use the pseudo-equivalent width of 
the Mg triplet lines at 5180 A as their third parame- 
ter, but we find that the reduced equivalent width of 
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Fig. 3. — (a) Color-magnitude diagram of observed stars in Segue 1. The large black circles represent the 71 stars identified as definite 
or probable radial velocity members of the galaxy using our subjective approach, the small black dots represent stars identified as likely 
or certain non-members, and the magenta crosses are spectroscopically confirmed background galaxies and quasars. The red curve shows 
the location of the red giant branch, subgiant branch, and main sequence turnoff populations in the globular cluster M92 and the cyan 
curv e shows the location of the horizontal branch of M13, both corrected for Galactic extinction and shifted to a distance of 23 kpc (data 
from lClem et aLl l2008l). (5) Spatial distribution o f obser ved stars in Segue 1. Symbols are the same as in (a), and the ellipse represents 
the half-light radius of Segue 1 from lMartin et al] 120081 ). (c) Velocity histogram of observed stars in Segue 1. Velocities are corrected to 
the heliocentric rest frame. The filled red histogram represents stars classified as members, and the hatched black-and-white histogram 
represents non-members. The velocity bins are 2 km s~ wide. 



the Ca triplet (CaT) lines works as well. iRutledge et alj 
(|1997al lbl) define this reduced equivalent width as W = 
ECa — 0.64(Vhb — V), where ECa is the weighted sum 
of the equivalent widths (EWs) of the CaT lines: SCa = 
0.5EW8498 + 1.0EW8542 + 0.6EW8662- Note that since 
the EM algorithm does not incorporate photometric in- 
formation, it must be run on a sample of stars that has 
already passed the color-magnitude selection. Also, EM 
may fail to select horizontal branch (HB) stars because 
their broad hydrogen lines interfere with measurements 
of the CaT EW; fortunately, the HB members of Segue 1 
are obvious. 

We display the distribution of observed stars in ra- 
dius, velocity, and reduced CaT EW in Figures S] and [5l 
Segue 1 stands out as the large overdensity of stars with 
velocities just above 200 km s~^ and smaller than average 
radii and W values. We caution that W is only a prop- 
erly calibrated metallicity indicator for stars on the red 
giant branch (RGB), which constitute a small minority 
of the data set examined here. Nevertheless, experiments 
with globular cluster stars reaching several magnitudes 
below the main seq uence turnoff from the compilation of 
iKirbv et al.l ()2010t ) show that while W does increase at 
constant metallicity toward fainter main sequence mag- 
nitudes, this increase is less than 2 A for stars within 2 
magnitudes of the turnoff. We therefore conclude that in- 
cluding both RGB and main sequence stars may broaden 
the Segue 1 distribution towards higher values of W 
(perhaps accounting for the clear presence of Segue 1 
stars in Figure [4)d up to W « 6 A), but should not 
significantly affect the performance of EM. 

The EM algorithm selects 68 stars as definite members 
of Segue 1 (membership probability p > 0.9). An addi- 
tional 3 stars have 0.8 < p < 0.9 and are classified as 
members by eye, yielding 71 very likely members. These 
71 stars correspond to 70 of the 71 subjectively classified 
membersjij Finally, 3 stars have membership probabili- 



The 



EM algorithm includes one star 



ties of 0.5 < p < 0.8. 

Our final results are based on a Bayesian analysis that 
allows for both contamination by Milky Way foreground 
stars and the contribution of binary orbital motions to 
the measured velociti es. These calculation s are a natural 
generalization of the iWalker et all (|2009b( ) EM method. 
The method is described in more detail in Paper II and 
is summarized here in §[3 In this framework, we find 
53 definite members {{p) > 0.9) and 9 further probable 
members (0.8 < (p) < 0.9), plus the 2 RR Lyrae variables 
(see S I4.2|) . but 7 of the stars considered likely members 
by the other two techniques receive lower probabilities 
of 0.4 < (p) < 0.8 here. With the exception of the dis- 
cussion in §[131 where we mention the range of velocity 
dispersions that can be obtained for different member 
samples, the main results of this paper (including the 
velocity dispersion, mass, and density of Segue 1) rely 
on this Bayesian analysis. It is important to note that 
unlike previous studies, we include all stars that pass 
our photometric cuts in the Bayesian calculations, not 
just the ones with high membership probabilities. Each 
star is weighted according to its probability of being a 
member of Segue 1 . This approach allows us to account 
correctly for the significant number of stars with mem- 
bership probabilities that are neither close to zero nor 
close to one. 

4. METALLICITY AND KINEMATICS OF SEGUE 1 

4.1. Stellar Metallicities and the Nature of Segue 1 

One of the defining differences between galaxies and 
globular clusters is that dwarf galaxies universally ex- 
hibit signs of internal chemical evolution and contain 

(SDSSJ100711. 80+160630.4) with a velocity of 247.1±15.9 km s-^ 
that we rated as too far removed from the systemic velocity to be a 
member candidate, and gives one star (SDSSJ100622. 85-1-155643.0) 
with a closer velocity but an even larger uncertainty 
(^hei = 223.9 it 37.8 km s^^) a lower membership probabil- 
ity of p = 0.73. Both of these stars are discussed further in the 
Appendix. 
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Fig. 4. — (a) Distribution of observed stars in velocit y an d radius. Filled red points represent stars that pass the color and magnitude 
selection (at either high or low priority) described in ii l2.2l and open black points are stars that lie outside that selection region. Stars 
that have been observed multiple times are plotted with their weighted average values. Segue 1 stands out as the large overdensity of stars 
near ti^gj = 200 km s~^ extending out to a radius of ^ 13'. Based on the distribution of Milky Way stars, it is clear that at small radii 
(r < 7') the level of contamination of the Segue 1 member sample is very low. In addition to Segue 1, there is also a distinct concentration 
of stars near 300 km s^^. (6) Distribution of observed stars in velocity and reduced Ca triplet equivalent width, a proxy for metallicity. 
As in the left panel, a large fraction of the Segue 1 members separate cleanly from the Milky Way foreground population. At W' > 5 A, 
the distributions begin to overlap, and unambiguously classifying individual stars as members or nonmembers becomes more difficult. 
Fortunately, relatively few stars are located in this region. It is clear that Segue 1 is more metal-poor than the bulk of the foreground 
population, although W' is a much less accurate metallicity indicator for main sequence stars than giants. The 300 km s^^ structure 
appears to be more enriched than Segue 1. 
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stars with a range of metallicities, while globulars gener- 



ally do not. Recent work has shown that multiple stellar 
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populations with different chemical abundance patterns 
are in fact present in some globular clusters, but these 
differences tend to be subtle (which is why they are only 
being recognized now) and are preferentially found in lu- 
minous clusters that are often argued to be the remnants 
of tidally stripped dwa rf galaxies (e.g.. iLee et alJ 119991: 
Marino ctal. 2009; Da Costa et~aLl 120091: iFerraro et ail 
2009; .Cohen et al..,20iaV 



We use the spectral synthesis metho d introduced by 
iKirbv. Guh athakurta. fc SnedenI (|2008l ) and refined by 
iKirbv et"ai. (2009, 20ld[) to measure iron abundancetQ 
in Segue 1 d i rectly from our medium-resolution spectra. 
IKirbv et alJ (|2010f ) showed that the metallicities mea- 
sured in this way are reliable for stars with log 5 < 3.6. 
Thus, we can only determine metallicities for the 6 red 
giant members of Segue 1; the fainter stars are all at or 
below the main sequence turnoff. The metallicities of 
these 6 stars span an enormous range, with two stars 
at [Fe/H] > -1.8 and two others at [Fe/H] < -3.3 
(see Table H]). One of the two extremely metal-poor 
(EMP) stars does not have a well-defined metallicity 
measurement because no Fe lines are detected in its spec- 
trum. The upper limit on its metallicity therefore de- 
pends on the assumptions, but it is certainly well below 
[Fe/H] = — 3. The mean metallicity of the Segue 1 red 
giants is [Fe/H] = —2.5, compa rable to the most m etal- 
poor galaxies identified so far (jKirbv et al.l 12008") . and 
the standard deviation, while not well-constrained with 
such a small sample, is ^ 0.8 dex. Using a completely 
independent data set, iNorris et al.l (|2010bD reach essen- 
tially identical conclusions regarding the Segue 1 abun- 
dance range and i dentify yet another EMP member star 
at [Fe/H] = -3.5 (INorris et al.l[20T0ah . 

The very large star-to-star spread in metallicities and 
the presence of EMP stars with [Fe/H] < —3.0 each 
independently argue that Segue 1 cann ot be a globular 
cluster contrary to initial s uggestions (jBelokurov et al.l 
I2007al: iNiederste-Ostholt et al. 2009). Only w Centauri 
among glob ular clusters ha s a comparable metallicity 
spread (e.g.. lNorris fc Da Cos ta 1995; Hilkcr ct al. 2004]), 
and that object is wi dely regarded to be the remnant 
of a dwa rf galaxv (iLee et al.l Il999t iMaiewski et al 



ICarraro fc Lial 120001: JHilker fc Richtlei 
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2OIOI ) rather than a true globular. 



The lowest metallicity Segue 1 giants are also at least a 
factor of four more met al-poor tha n any known star in 
a glo b ular cluster (e . g.. IKing et al.llT9 98: Kraft fc IvansI 
[200I IPreston et all I2OO6I: ICarretta et al.l l2009l ). We 
conclude that the metallicities of stars in Segue 1 pro- 
vide compelling evidence that, irrespective of its current 
dynamical state. Segue 1 was once a dwarf galaxy. 

4.2. Binary and Variable Stars 

Before attempting to determine the velocity dispersion 
and mass of Segue 1 , we consider the impact of variable 
stars and binaries that could be in our sample. The G09 

^^ For historical reasons, these c alculations use a solar iron abun- 
dance of 12 + loge(Fe) = 7.52 IIAnders fc Grev cssc 1989), but 
the difference between this assumption and the modern value of 
12 -I- log e(Fe) = 7.50 ± 0.04 HAsplund et al.|[2009D is negligible. 



members include two horizontal branch stars in Segue 1. 
Our repeated measurements demonstrate that the ve- 
locities of both of these stars vary with time, leading 
us to conclude that they are RR Lyrae variables. Fol- 
lowup photometry with the Pomona College 1 m tele- 
scope at Table Mountain Observatory confirms that one 
of these stars, SDSSJ100644.58-H155953.9, is a photo- 
metric variable with a characteristic RR Lyrae period of 
0.50 daysQ We do not detect variability in the second 
star, SDSSJ100705.60-I-160422.0, but the limits we can 
place are not inconsistent with the low amplitude vari- 
ability that might be expected for such a blue star. Given 
their blue colors, the stars are probably type c variables 
pulsating in the first overtone mode. Because the light 
curve phases at the times our spectra were acquired are 
not known, we cannot measure the center-of-mass veloc- 
ities of these stars and must remove them from our kine- 
matic sample even though they are certainly members of 
Segue 1. 

We also obtained multiple measurements of five of the 
six Segue 1 red giants in order to check whether any of 
them are in binary systems. For four of the stars, the 
velocity measurements do not deviate by more than 2 a 
from each other, providing no significant indication of bi- 
narity, although long period or low amplitude orbits can- 
not be ruled out. SDSSJ100652.33-I-160235.8, however, 
shows clear radial velocity variability, with the velocity 
decreasing from 216.1 ± 2.9 km s~^ on 2007 November 
12 to 203.0 ± 2.3 km s^^ on 2009 February 27, and then 
rising back to 210.8 ± 2.3 km s^^ on 2010 February 13. 
Assuming that these observations correspond to a single 
orbital cycle, we infer a period of ~ 1 yr and a companion 
mass of - 0.65 + 0.25(l/sin3i ~ 1) Mq. 

With at least one out of the brightest six stars (ex- 
cluding the even more evolved horizontal branch stars) 
in the galaxy in a binary system, the binary fraction 
of Segue 1 is likely to be s ignifica nt, as has b een found 
for other dwarf galaxies (IQueloz . Dubath, fc PasauTml 
119951 : lOlszewski. Prvor. fc Armandrofil [1996) and some, 
altho u gh not all, globular clusters (Fischer et al.l 
19931 JYan fc CohenI 119 96: Rub enstein fc Bailvnl 11997 



Clark. Sandauist. fc Boha .2004 iSolhma et al.l [20071) . 
For main sequence stars, which dominate our sam- 
ple, the binary fraction can only be larger since 
the tight binary systems will not yet have been de- 
stroyed by the evolution of the more massive com- 
ponent. We have a limited sample of repeat obser- 
vations of some of the main sequence members, in 
which two additional stars, SDSSJ100716. 26-^160340.3 
and SDSSJ100703.15-fl60335.0, are detected as proba- 
ble binaries. However, these binary determinations are 
almost certainly quite incomplete, and proper corrections 
for the inflation of the observed velocity dispersion of 
Segue 1 by binaries must be done in a statistical sense, 
as we discuss in §[5j 

4.3. Kinematics 

Because the issues of membership and binary stars are 
so critical to our results, we must experiment with dif- 

^^ Periods of 0.5 d are of course subject to the possibility of 
aliasing, and the light curve is not complete enough to rule out a 
period of 1.0 d. However, such long periods are extremely rare for 
RR Lyraes, so we consider the 0.50 d period to be the most likely 
solution. 
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ferent samples of member stars and methods of deter- 
mining the velocity dispersion. Inspection of Figure |4] 
makes clear that for W < 3 A, the expected contamina- 
tion by Milky Way foreground stars is negligible (1 — 2 
stars). Very conservatively, then, we can select the stars 
with W^' < 3 A and 190 km s'^ < u < 225 km s"^ as 
an essentially clean member sample (the exact velocity 
limits chosen do not matter, since the next closest stars 
are at w = 175 km s~^ and v > 300 km s~^). With the 
two RR Lyrae variables and the one obvious RGB binary 
removed, the velocity dispersion of the other 34 stars is 
3.3 ± 1.2 km s~^. (Note that we calculate the velocity 
dis persion using a maximum likelihood method follow- 
ing I^MSl^ES! 12006.) This value can be regarded in 
some sense as a lower limit to the observed dispersion of 
Segue 1 (prior to any correction for undetected binaries) . 

Since this conservative approach involves discarding 
nearly half of the data, we would also like to consider 
alternatives. The largest member sample that we can de- 
fine is the 71 stars selected using either our holistic, sub- 
jective criteria in ? l3.1l or the EM algorithm. After again 
excluding the two RR Lyraes and the RGB binary, the 
raw velocity dispersion of these stars is 5.5 ± 0.8 km s~^, 
which we take as an upper limit to the observed disper- 
sion of Segue 1 (as before, prior to correcting for unde- 
tected binaries). However, the dispersion in this case is 
dominated by a single star (SDSSJ100704.35+160459.4; 
see § l4.4p . If we remove this object from the sample, the 
dispersion of the remaining stars falls to 3.9 ±0.8 km s^^, 
corresponding to a factor of 2 decrease in the mass of the 
galaxy. Other stars that are borderline members or non- 
members have negligible effects on the derived dispersion 
(see Appendix). 

We note at this point a curious finding regarding the 
brightest stars in Segue 1. If we isolate the evolved stars 
(6 giants and 2 horizontal branch stars) in the sample, 
their velocity dispersion appears to be quite small. For 
the HB stars and the binary on the giant branch we can- 
not assume that we have enough measurements to aver- 
age out the effects of the binary orbit and RR Lyrae pul- 
sations, but we estimate a dispersion of 1.31q7 ^™ ^^^ 
for the other five RGB stars. Using our full Bayesian 
analysis (§ [5]) and including the binary, the intrinsic dis- 
persion of the giants is 2.0J^^'7 km s^^. Given the sub- 
stantial error bars, these values are formally consistent 
with the larger dispersion obtained for the full data set, 
even though they are also close to zero. Nevertheless, 
the velocity dispersion we determine for the remaining 
stars is not significantly affected by the inclusion or ex- 
clusion of the giants and horizontal branch stars. With- 
out any known physical mechanism that could change the 
kinematics of Segue I for stars in different evolutionary 
states, we conclude that the apparently small dispersion 
of these stars is most likely a coincidence resulting from 
small number statistics. 

Next , we use the sample defined by the iWalker et al.l 
()2009b[) EM algorithm. Since this sample is nearly iden- 
tical to that considered in the previous paragraphs, the 
results are unchanged: a dispersion of 5.7±0.8 km s~^ for 
all 71 stars minus the RR Lyraes and the RGB binary, 
and 4.1 ± 0.9 km s"^ when SDSSJ100704.35-I-160459.4 
is removed. It is worth noting that all of these measure- 
ments and those described above are consistent within la 



with the original velocity dispersion of 4.3 ± 1.2 km s~^ 
determined by G09. 

The two above methods make various assumptions 
about how membership is defined, but do not allow for 
a fully self-consistent statis tical treatment. The method 
described in Paper II (iMartinez et al.|[2010l ) and summa- 
rized in §[5] treats these assumptions and the data analy- 
sis in a fully Bayesian manner. This analysis identifies a 
total of 61 stars (excluding the 2 RR Lyrae variables) as 
likely members with (p) > 0.8. This method arrives at 
significantly lower membership probabilities for 8 stars 
compared to the EM and subjective analyses. These 
stars fall into three partially overlapping categories: ve- 
locity outliers, frequently with large velocity uncertain- 
ties as well (such that there is a significant chance that 
the star's true velocity is far away from the systemic ve- 
locity of Segue 1); stars with large reduced CaT EWs 
{W > 4 A); and stars at large radii (r > 10'). The only 
one of these stars that has an appreciable effect on the ve- 
locity dispersion is SDSSJ100704.35-I-160459.4. Remov- 
ing each of the 7 other stars from the sample changes 
the dispersion by less than 0.2 km s~^. With only the 61 
most likely members included in the calculation, we find 
a velocity dispersion of 3.4 ± 0.9 km s^^. 

4.4. Individual Stars With Ambiguous Membership 

Despite our best efforts to define a member sample that 
is both clean and complete, there are fundamental uncer- 
tainties that cannot be avoided in determining whether 
any given star is a member of Segue 1. In particular, we 
know that Segue 1 is an old, metal-poor stellar system lo- 
cated 23 kpc from the Sun, and moving at a heliocentric 
velocity of ^ 207 km s^^. Unfortunately, the Milky Way 
halo also contains old, metal-poor stars that span ranges 
in distance and velocity that encompass Segue 1. Given 
a large enough search volume, it is therefore inevitable 
that some halo stars with the same age, metallicity, dis- 
tance (and hence the same colors and magnitudes), and 
velocity as Segue 1 will be found. We can use observa- 
tions and models to estimate the expected number of such 
stars included in our survey, but that does not help us in 
ascertaining the provenance of individual stars. 

As a result, we are left with a small number of stars 
whose membership status is necessarily uncertain. The 
algorithms discussed in § 13. H allow us to assign member- 
ship probabilities to these objects, which is the best sta- 
tistical way to deal with our limited knowledge. Never- 
theless, each star that we observed either is or is not a 
member, and with a small sample, the assumption that 
a particular star has a membership probability of, e.g., 
0.6 can produce different results than if it were known 
absolutely to be a member or not. Most of the stars in 
this category do not have an appreciable effect on the 
derived properties of Segue 1 (most importantly the ve- 
locity dispersion), either because they lie near the middle 
of the distribution or because their velocity uncertainties 
are relatively large. One object, however, can make a 
significant difference, as discussed in the following para- 
graph. A few other stars that cannot be classified very 
confidently in one category or the other are listed in the 
Appendix, but their inclusion or exclusion has minimal 
impact on the properties of Segue 1 or the results of this 
paper. 
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The one star that can individually affect the kinemat- 
ics of Segue 1 is SDSSJ100704.35+160459.4. This star 
is a major outlier in velocity, with a mean velocity from 
two measurements of 231.6 ± 3.0 km s""'^. It is therefore 
located more than 6a (where a is estimated from the rest 
of the stars) away from the systemic velocity of Segue 1, 
but because its position is extremely close (38") to the 
center of the galaxy and it has a CaT EW that could 
plausibly be associated with Segue 1 (although on the 
high side), the EM method returns a membership prob- 
ability of 1. The membership probability from our full 
Bayesian analysis (see Paper II) is lower, but far from 
negligible, at 0.49. The high velocity of this star relative 
to the systemic velocity of Segue 1 could be explained if 
it is a member of a binary system, but our two velocity 
measurements of it (separated by 1 yr) do not show a 
significant change in velocity, so the period would have 
to be > 5 yr. In addition to its disproportionate effect on 
the velocity dispersion, for an equilibrium model a star 
that is a 6a outlier from the mean velocity but is located 
so close to the center of the galaxy implies strongly radial 
orbits. None of the other stars in the sample lead to a 
preference for extreme velocity anisotropy. 

5. THE INTRINSIC VELOCITY DISPERSION AND 
MASS OF SEGUE 1 

Having carefully considered the membership of each 
star and the effects of the key outliers in the previous two 
sections, we are now in position to determine the best 
estimate of the intrinsic velocity dispersion of Segue 1 
based on all available data. The full c alculations that 
we use for this purpose are presented in Martin ez et al.l 
(|201Q) . but we include a summary here for convenience. 

5.1. Binary Correction Method 

Given a sample of stars that may be members of 
Segue 1, we allow for the possibility that some of the 
observed stars are likely members of binary star systems 
rather than single stars. We therefore must treat the (un- 
known) velocity of the star system's center of mass Vcm 
and the measured velocities themselves as distinct quan- 
tities. For each star system of absolute magnitude My 
we have a set of measured velocities Vi, where i runs over 
the number of repeat observations of the star under con- 
sideration at times ti, and the associated measurement 
uncertainties e^. We write the likelihood of obtaining the 
observed data for each star assuming it is a member of 
Segue 1 (SI) in terms of a joint probability distribution 
in the measured velocities Vi and the unknown center-of- 
mass velocity Vcm- 

/:sl(^'^|e^,i^,My;CT,^,B,^) 

P{v^\Vcm,et,ti, My; B, ^)P{Vcm\<J, l^)dVcm{^) 



cx(l-B) 



(M-fil^ 



V^^+^ 



+ BJ{a,fi\^), 



where the first factor in the integrand is the probability 
of drawing a set of velocities Vi given a center-of-mass 
velocity Vcm and certain values for the binary param- 
eters B and 3^. B represents the binary fraction and 



^ is the set of binary parameters [/iiogPjCriogp] (see be- 
low) . The second factor is the probability distribution of 
center-of-mass velocities given an intrinsic velocity dis- 
persion of Segue 1, a, and a systemic velocity, /i. In 
the last line of Equation [l] (v) is the average velocity 
weighted by measurement errors and am is the uncer- 
tainty on this weighted average for the combined mea- 
surements of each star. Note that we assume that the 
center-of-mass velocity distribution of Segue 1 is Gaus- 
sian. We also use metallicity (W) and position to help 
determine membership, so that the full likelihood is of 
the form C{vi,W' ,r), but we omit the metallicity and 
position dependence in the equations presented here for 
simplicity. The absolute magnitude of each star is taken 
into account so that its radius can be calculated and 
only binaries with separations l arger than the stella r ra- 
dius are allowed, as described in lMinor et al.l (|2010( ) and 
Paper II. 

For each star, the J(/i, a|^) factor is generated by 
running a Monte Carlo simulation over the distribution 
of binary properties, which include the orbital period, 
mass ratio, and orbital eccentricity. Unfortunately, the 
characteristics of binary populations in dwarf galaxies 
are completely unknown at present. The best empirical 
constraints on binary properties come from studies of the 
Milky Way, but there is no guarantee that the small-scale 
star forming conditions in the Milky Way and those that 
prevailed in Segue 1 many Gyr ago are similar, especially 
since Segue 1 has a metallicity two orders of magnitude 
below that of the Galactic disk population. In princi- 
ple, a lower metallicity could change the probability of 
forming binary systems (and higher order multiples) and 
the properties of those systems, but binaries with sep- 
arations in the range that can affect our observations 
(< 1 AU) are ex pected to form via disk fragmentation 
(e.g.. lKratter et al. 2010) or interactions between proto- 
stellar cores (e.g.. .BatG.2009l ). neither of which should be 
very sensitive to metallicity (M. Krumholz 2010, private 
communication) . 

In the Monte Carlo simulations, we fix the distribu- 
tions of the mass ratio and eccentricity to follow that 
observed in solar neighborhood binaries. However, since 
the period distributions of binary populations have been 
obser ved to differ dramatically f rom cluster to cluster 
(e.g., iBrandner fc Koehleii 119981: iPatience et all [2002), 
we allow for a range of period distributions. Specifi- 
cally, we assume the distribution of periods has a log- 
normal form similar to that of Milky Way field bin aries 
puauennov fc Mavoi|[l99lHRaghavan et al. II2010D . We 
take the mean log period /xiogp and the width of the 
period distribution criog p as free model parameters. We 
further assume that the period distribution observed in 
Milky Way field binaries is the result of superposing nar- 
rower binary distributions from a variety of star-forming 
environments. Thus, our prior on the period distribu- 
tion is that it is nar rower than, but consistent w ith be- 
ing drawn from, the iDuquennov fc MavoH (|1991h period 
distribution with /iiogp — 2.23 and aiogp = 2.3. To ac- 
complish this, we choose a fiat prior in criog p over the 
interval [0.5, 2.3] and a Gaussian prior in /iiogp centered 
at /iiogp = 2.23, with a width chosen such that when 
a large number of period distributions are drawn from 
these priors, they combine to reproduce the Milky Way 
period distribution of field binaries. We can then write 



Complete Spectroscopic Survey of Segue 1 



11 



the likelihood of each star being a member of Segue 1 or 
the Milky Way as 



£(?;,|e„i,,My;/,B,(T,/x,^) = 
(1 - f)CMw{v^\e^) + fCsi{v.,\e,, t,,Mv; B, a, /x, ^p) 

where / is the fraction of the total sample that are 
Segue 1 members. The second term in Equation [5] is 
given by Equation [1] and the likelihood of membership 
in the Milky Way is 



CMwivi\ei) ex 



^^^ 



V2^ 



-PMw{Vcm)dVcm- (3) 



The expected Milky Way velocity distributio n is taken 
from a Besangon model (|Robin et alJ l2003| ) after the 
same photometric criteria used to select Segue 1 stars 
have been applied. 

5.2. Binary Correction Results 

By applying the above analysis to our full Segue 1 data 
set we can correct for the likely presence of binaries in 
the sample and derive the intrinsic velocity dispersion 
of the galaxy. All observed stars (not only the likely 
members) that meet the photometric selection cut de- 
scribed in ii l2.2l are included in this calculation, with the 
weight for each star determined by its likelihood of mem- 
bership in Segue 1. Our posterior probability distribu- 
tion for the dispersion is maximized at ct = 3.7 km s~^, 
~ 12% smaller than what we measure without a binary 
correction. The la uncertainties on the dispersion are 
-1-1.4 km s^^ and —1.1 km s~^. We find a 90% lower limit 
on the dispersion of 1.8 km s~^, and the probability of the 
true velocity dispersion being less than 1 km s~^ is ~ 4%. 
The differential and cumulative probability distributions 
for a are displayed in Figure |6l If the gravitational po- 
tential of Segue 1 were provided only by its stars, the 
velocity dispersion would be < 0.4 km s~^ (G09). Our 
lower limit on the dispersion therefore allows us to con- 
clude with high confidence that Segue 1 is dynamically 
dominated by dark matter unless it is currently far from 
dynamical equilibrium. We discuss the implausibility of 
large deviations from equilibrium caused by tidal forces 
later in §[Ll 

We further note that small intrinsic velocity disper- 
sions (cr < 2 km s~^) can only be obtained if the binary 
period distribution is skewed toward short periods. In 
particular, mean periods of less than ~ 40 yr are re- 
quired to produce such a small dispersion (Paper II) ; for 
comparison, t he mean period in the solar neighborhood 
is 180 yr (Duquennov fc Mavorill99l"l ). The posterior dis- 
tribution we derive for the mean period in Segue 1 is in 
fact weighted toward quite short periods (~ 10 yr). This 
result is not an artifact of the short time baseline of our 
multi-epoch data (observations on timescales of ~ 1 yr 
cannot constrain periods of centuries or longer), because 
our priors allow for flatter period distributions than the 
posterior for the mean period. Despite this preference for 
shorter periods, the data strongly indicate only a modest 
contribution to the velocity dispersion from binaries be- 
cause the period distribution is still wide and the small 



number of detected binaries indicates that the fraction 
of stars in close binary systems is not large. 

To provide reassurance that our priors on the period 
distribution are not biasing the posterior period distri- 
bution toward longer periods (and hence the corrected 
velocity dispersion toward higher values), we repeated 
the calculations above with a flat prior on the mean pe- 
riod /xiogp. The best- fit mean period barely changed, 
and while somewhat shorter periods are allowed in this 
case, the probability of an intrinsic velocity dispersion 
less than 1 km s^^ does not increase (see Figure[7]). Even 
with a more extreme logarithmic /iiogp prior, the likeli- 
hood of a small dispersion is unchanged despite the re- 
sulting very short mean period. The reason for this out- 
come is that when the mean period is forced to be short, 
the binary fraction is then constrained to be low and the 
width of the period distribution is similarly constrained 
to be small to fit the observed changes in the velocities 
and the observed velocity distribution. Hence, the tail 
of the probability distribution toward low velocity dis- 
persions cannot be made significantly larger by having a 
prior that biases the result to shorter mean periods. 

These results co nt rast with the findings of 
IMcConnachie fc Cotel (|2010t ). who conclude that 
galaxies like Segue 1 could have very low intrinsic 
velocity dispersions that are inflated substantially by 
the presence of binary stars. While we agree with their 
results given the assumptions they make, two primary 
factors are responsible f or the different conc lusions 
from our analysis. First, IMcConnachie fc Cotel ignore 
binaries with periods longer than 10—100 yr. This cutoff 
appears reasonable from an observational perspective, 
since such binaries will not be detectable in current 
data sets, but it has the effect of making the binary 
fractions they require very large because the majority 
of Milky Way binaries have periods longer than 100 yr. 
Second, they analyze only single-epoch velocity data 
sets, whereas the multiple measurements we have for a 
number of stars give us much greater leverage with which 
to determine the inflation of the velocity dispersion 
caused by binaries. Our Bayesian analysis including 
the multi-epoch data and all the information in the tail 
of the velocity distribution shows that a substantial 
inflation by binaries is disfavored for the Segue 1 data 
set presented here (see Paper II for more details). 

5.3. Mass of Segue 1 

The same Bayesian machinery described in i) l5.1l for de- 
termining the velocity dispersion of Segue 1 can also be 
employed to calculate the mass of the galaxy (again in- 
cluding a corr ection fo r binary stars) . IWolf et aLl (J20101 
also see Walk er et al.l [2009a) derived a simple formula 
for the mass within the half-light radius of a system: 
Ml/2 = 3(7^ri/2/G. For a fiat prior on a (see Paper II 
for a discussion of the effect of the choice of priors), 
we find a posterior probability distribution on the mass 
within the 3D half-light radius of Segue 1 (38 pc) of 
5.8^3 ^ X 10^ M0, consistent with the mass determined 
by G09 from the original data setEj Because the uncer- 
tainties on the mass are not Gaussian, this measurement 

^"^ Note that calculating A/j ;2 directly from the stellar velocity 
data set with this Bayesian approach is not the same as sim ply 
plugging the derived values of a and r^M into the IWolf et al.l for- 



12 



Simon et al. 



c 
-a 



o 



1.0 






/\ 


\ 


. , 













/ / 


\ 


\ 


binary- 


- 








/ / 


\ 


\ 


corrected 


- 


0.8 


- 




/ / 




\ \ 


no binary 


- 










correction 


- 








/ / 




\ \ 




; 


0.6 


- 




/ / 




\ \ 




- 








/ / 




\ \ 




. 


0.4 


- 


1 


/ / 
/ / 
/ 




\ \ 
\ \ 
\ \ 




- 


0.2 


- 


/ 


1 




\ \ 




- 






/ 




^^ 


\ 


- 






^/ 


/ 




\ 


s. N 


- 


on 


■^ 


< 






1 


^^..^^^ 





2 4 6 8 

Velocity dispersion (km s~') 



10 



^ 1.0 


1 . , , 1 . , , 1 ^^„^i^ 1 


b 


y^"^ 


V 


y^ 


^— ' 


/ 


>, 0.8 


- / - 


+j 
















XI 


/ 


« 0.6 


1 


o 


j 


!-, 


/ 


a. 


/ 


a, 0.4 


- 1 


> 


/ 






.1-) 




ffl 


/ 


3 0.2 


/ 


e 


/ 


3 


y^ 


" 0.0 


— — -^^ . , , 1 . , , 1 , . , 1 , , . 



2 4 6 8 

Velocity dispersion (km s~') 



10 



Fig. 6. — (a) Posterior probability distribution for Segue 1 velocity dispersion before (dashed) and after (solid) correcting for binary 
stars. (6) Cumulative probability distribution for Segue 1 velocity dispersion after correcting for binaries. 
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binary period on the derived velocity dispersion. The black 
solid curve shows on the Segue 1 velocity dispersion for our 
preferred assumption of th e Milky Way field binary prior from 
IDuouennov & Mavoii 11199 II ). The dashed red curve represents a 
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the result of a logarithmic prior that is even more strongly biased 
toward short periods. The very small changes in both the most 
likely value of the velocity dispersion and the size of the tail to 
low values of the dispersion (ct < 1 km s~^) demonstrate that 
our results are robust to differing assumptions about the binary 
population in Segue 1. 



disagrees with the stellar mass of Segue 1 (^ 1000 M©) 
at much more than 1.8ct significance; the 99% confidence 
lower limit on the mass is 16000 M© (however, the mag- 
nitude of this low-mass/low-cr tail in the probability dis- 
tribution is prior-dominated) . The V-band mass-to-light 
ratio within the half-light radius is ^ 3400 Mq/L©. Since 
there is no evidence that the dark matter halo of Segue 1 
is truncated at such a small radius, this value represents a 
lower limit on the total mass-to-light ratio of the galaxy, 
which could be 1-2 orders of magnitude larger. 

mula. The final value for M^n is very similar, but the uncertain- 
ties are more accurately determined (in particular the 2 and Scr 
confidence intervals) when we have determined the full probability 
distribution. 



6. A STREAM AT 300 km s^^ 

In addition to Segue 1 and the Milky Way foreground, 
we clearly detect a third population of stars in our kine- 
matic data, at a heliocentric velocity of 300 km s~^. We 
refer to this structure as the "300 km s~^ stream" be- 
cause of its lack of spatial concentration within our sur- 
vey area. However, we recognize that these stars could 
still be part of a bound system as long as the angular size 
of the object is comparable to or larger than our field of 
view (diameter ^ 20', which corresponds to a physical 
size of at least 116 [(i/20 kpc] pc). As an example, a 
galaxy similar to And XIX, which has a half-light ra- 
dius of ^ 1.7 kpc, would subtend s everal degrees at this 
distance (jMcConnachie et al.ir2008f) . 

G09 also recognized the existence of this stream, find- 
ing 4 stars in it among their smaller sample. We now 
present conclusive confirmation that this structure is 
real, with ~ 20 stars in our new data set (see Figures H]- 
[5]). Because of the smaller member sample and the low 
contamination from Milky Way stars at such extreme 
velocities, the probabilistic membership algorithms de- 
scribed in §[3TT] are not necessary in this case. Instead, 
we select stars that have velocities between 275 km s~^ 
and 325 km s~^ (the exact velocity limits are not im- 
portant; see Fig. |8}:) and meet the same photometric cri- 
teria that were used for Segue 1. The color/magnitude 
screen eliminates 5 stars, leaving 24 likely members in 
the stream. 

The CMD of these stars is similar to that of Segue 1, 
indicating that the stream is also ~ 20 kpc away (see Fig- 
ure [8]). The stream main sequence appears to be slightly 
redder (suggesting a higher metallicity) and slightly 
closer than Segue 1, although any differences are near 
the limit of what can be determined from the SDSS data. 
The stream stars are matched quite well with the fiducial 
sequ ence of th e glob ular cluster M 5 ([Fe/H] — —1.27) 
from lAn et al.l ()2008() . supporting the higher metalHcity 
that one would have guessed by eye. By comparing to 
various globular cluster sequences, we estimate a distance 
of sa 22 kpc and a metallicity of [Fe/H] « —1.3, but the 
uncertainties on both numbers are substantial. 

The 24 candidate members have a mean veloc- 
ity of 298.8 ± 1.7 km s~^ and a velocity disper- 
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Fig. 8. — (a) Color-magnitude diagram of observed stars in the 300 km s~^ stream. The filled black circles represent stars identified as 
candidate stream members, while the open gray circles are the Segue 1 members. (6) Spatial distribution of observed stars in the stream. 
Symbols are the same as in (a), (c) Velocity histogram of observed stars in the stream. Velocities are corrected to the heliocentric rest 
frame, and the velocity bins are 2 km s~^ wide. 



sion of 7.0 ± 1.4 km s~^, comparable to the disper- 
sion of other known stream s (Chap man et al.l 120081: 
Grillmair. Carlin. fc Ma icw skil [20081 : lOdenkirchen et ahl 
2009; Ncwbcrg et al. 2010). Several of these stars may be 
foreground contaminants; in particular, the bright star 
at g-i = 0.48, r = 17.60 (SDSSJ100720.00-I-160137.5) 
is located a bit below the horizontal branch if the 
distance of 22 kpc preferred by the main-sequence 
fitting is used, although it does still lie within the 
AGB selection region shown in Figure [1] However, if 
the stream is at a slightly larger distance then this 
star could well be a horizontal branch member. The 
bluest of the faint stars (SDSSJ100650.83-hl60351.2; 
g — i = 0.28, r — 21.97) is ~ 3a away from the 
fiducial sequence used for the original target selection, 
and two other stars (SDSSJ100732.48-I-160500.5 and 
SDSSJ100708. 38-1-155646.3) have reduced CaT EWs that 
are well above those of the bulk of the stream popula- 
tion. Even if we remove all four of these stars, the stream 
properties do not change significantly; the mean velocity 
in that case is 298.7 ± 1.5 km s^^ and the dispersion is 
5.6 ±1.2 kms-i. 

If we assume that the 300 km s^^ stream is a bound 
structure wit h a half-light rad ius larger than our sur- 
vey area, the iWolf et all ()2010[ ) formula implies a lower 
limit on the mass contained within its half-light radius 
of 5.3 X 10^ (r/116 pc) Mg. This value would place the 
stream on the mass-radius relation of Milky Way dwarf s 
for a 3D half-light radius of ~ 500 pc (IWolf et al.ll2'010h . 
On the other hand, it may be worth noting that the 
stream appears to be more extended in the east-west di- 
rection than north-south (cf. Figures [3^ and [8^), consis- 
tent with an east-west extent. We therefore tentatively 
suggest that the 300 km s^^ stream could be the kine- 
matic counterpart of the similarly oriented p h otom etric 
feature identified by iNiederste-Ostholt et al.l (|2009l see 
^ 17.1. ip . An alternative possibility is that the stream 
might be related to Leo I, which is located approxi- 
mately 3.8° due sou th of Segue 1 at a similar velocity 
(282.9 ± 0.5 km s'^ iMateo. 01sz"CT^ki. fc Walker 2008 .) 
and metalhcity ([Fe/H] = -1.20: iKirbv et al.l l20ill ). 
Tracing the stream with wider field spectroscopic data 
to test these hypotheses would be very desirable. 



7. IS SECUE 1 UNDERGOING TIDAL 
DISRUPTION? 

INiederste-Ostholt et al.l (|2009l hereafter NO09) argued 
that rather than being a bound, dark matter-dominated 
dwarf galaxy. Segue 1 is a tiny star cluster whose appar- 
ent velocity dispersion has been inflated by contamina- 
tion from the Sagittarius (Sgr) stream. We have demon- 
strated that Segue 1 is not a star cluster (ii l4.1|) . but that 
flnding does not address the issues of tidal disruption or 
contamination. However, our new observations and a re- 
analysis of the SDSS data present some difficulties for the 
NO09 hypothesis. NO09 base their argument on several 
key points: (1) over a large area around Segue 1, there are 
very low surface density features whose color-magnitude 
diagrams are very similar to that of Segue 1 itself; (2) the 
surface brightness profile of Segue 1 appears to depart 
from a standard King or Plum mer model at large radii; 
(3) the iFellhauer et al.l (|2006t ) model of the Sgr stream 
predicts that there should be some very old Sgr material 
near the position and velocity of Segue 1; (4) using SDSS 
blue horizontal branch (BHB) stars, a coherent feature 
that can probably be identified with the Sgr stream ap- 
proaches the position and velocity of Segue 1; (5) given 
the larger velocity dispersion of Sgr, only a small amount 
of contamination of the Segue 1 member sample by Sgr 
stars is necessary to substantially inflate the apparent 
velocity dispersion. 

Below we discuss each of these ideas in turn and con- 
sider how our results affect their interpretation. We show 
that: (1) while photometric tidal features are indeed 
present in this held, they appear more likely to be as- 
sociated with other kinematically detected tidal struc- 
tures rather than Segue 1; (2) in the radial range where 
the photometric and kinematic constraints are good, the 
surface brightness proflle is well described by a Plum- 
mer model; (3) observational evidence for older wraps of 
the Sgr stream is non-existent, and even if present, more 
recent models suggest that this material has a very low 
surface density and differs in velocity from Segue 1; (4) 
the BHB feature identifled by NO09 as potentially con- 
taminating the Segue 1 data set is offset noticeably in 
both position and velocity from Segue 1; (5) the result- 
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ing contamination has therefore been overestimated, and 
furthermore, these BHB stars seem to be associated with 
the Orphan Stream (in which case they are too spatially 
confined to affect observations near Segue 1) rather than 
the Sgr stream. We therefore conclude that contamina- 
tion by Sgr stream stars does not have a significant im- 
pact on the measured velocity dispersion of Segue 1. We 
then consider the evidence that Segue 1 could be tidally 
disrupting, finding that while it is not possible to rule 
out recent tidal disturbances, the existing data do not 
provide significant support for such an interaction. 

7.1. Reconsidering the Disrupting Cluster Scenario 
7.1.1. Extended Tidal Debris Near Segue 1 

As NO09 have shown, there is no doubt that there are 
spatially extended structures whose stars roughly follow 
the Segue 1 fiducial sequence distributed over a wide area 
around Segue 1. Within ^ 1° of Segue 1, this population 
is even visible by eye in SDSS color-magnitude diagrams 
(CMDs). 

What is less obvious is that these stars are necessarily 
associated with Segue 1. One alternative is that they are 
instead part of the Sgr stream. After all, it is clear from 
both observations (Bclokurov ct aL 2007a) and models 
(G09, iLaw. Johnston, fc Maiewskil [2005] ) that the Sgr 
stream passes through this part of the sky at a dis- 
tance similar to that of Segue 1. Indeed, at the po- 
sition of Segue 1, th e stream runs nearly east-west 
(jBelokurov et al.ll2006( ) , exactly matching the orientation 
of the features identified by NO09. Another possibility is 
that the tidal features could be related to the 300 km s"^ 
stream, which also shares a very similar stellar popula- 
tion to Segue 1 (note that at the relevant distances, most 
of the stars detected by Sloan are on the main sequence, 
and thus the CMD filtering is primarily sensitive only 
to distance, not to metallicity). In either case, it seems 
more natural to associate this apparent tidal debris with 
one of the two known tidal structures at this position, 
rather than with the one object that is not obviously 
undergoing tidal disruption. 

7.1.2. The Surface Brightness Profile of Segue 1 

A second facet of the NO09 picture is the apparent 
excess of stars above the fitted Plummer, King, and 
exponential models at large radii. However, NO09 
themselves note that the area covered by their deeper 
imaging is not large enough to define a meaningful back- 
ground level, calling into question the significance of 
this excess. Within the radius probed by our kinematic 
data (~ 3 half-light radii), their photometric analysis 
shows that the data are fit well by a Plummer model, 
in agreement with the distribution of spectroscopically 
confirmed member stars that we derive (see § 17. 2|) . 
"Extratidal" excesses similar to the one claimed in 
the outer parts of Segu e 1 have b een seen in many 
other dwarf sp heroidals (llrwin fc H atzidimitriou 119951: 
Maiewski et al.l l2000aL 120051 : iMartincz-Dclgado et al ! 
20011 : iPalma et all 120031 : Walchcr ct al. 2003: 



Wilkinson ct a ll 120041: [M ashchcnko. Couchman. & Sills 



2005| :.WcstfaU et al.112006;: .Mufioz et al...200&: . Sohii et al.l 



7.1.3. Ancient Wraps of the Sagittarius Stream 

The Sagittarius dwarf has both leading and trailing 
tidal streams stretching across the entire sky. The most 
recent wrap of the strea ms has been detected robustly in 
numerous ways (e.g.. Ilbata et al.l 120011: I Vivas et aF 



2001 



iDohm-Palmer et al.l 120011: Bellazzini et al 



2003t iNewberg et all I200I iMaiewski et al.l 120031 : 



Belokurov et al.ll2006[ ). Models predict that Sgr debris 
stripped on previous orbits may be present as well, but 
there is currently little observational evidence for such 
material. G09 showed that while Segue 1 is spatially 
coincident with the leading arm of the Sgr stream, 
the velocities of the recently stripped stars differ from 
that of Segue 1 by ~ 100 km s~^, firml y ruling out an 
association. NO09 pointed out that the iFellhauer et all 
(|2006.) model predicts that there are also Sgr stars 
stripped several orbits earlier at this position that have 
vel ocities similar to Segue 1 . In the most recent model 
by ILaw fc Maiewskil (|2010al ) , which is the most success- 
ful to date in matching observations [3 however, stars in 
this ancient wrap uniformly have much lower velocities 
in this part of the sky {vqsr < 6 km s^^, compared 
to vgsr = 113 km s~^ for Segue 1). In addition, the 
surface density associated with the stars stripped at the 
earliest times should be overwhelmingly smaller than 
tha t of the more recent debr is. Within 10° of Segue 1, 
the ILaw fc Maiewskil ()2010aD simulation contains > 200 
times as many stars in the recently-stripped, leading 
(negative velocity) stream as are present in the older, 
trailing stream. Since NO09 find similar numbers of 
stars in their observed positive and negative velocity 
BHB streams, this provides a strong argument that 
the BHB structure closer in velocity to Segue 1 i s not, 
in fact, related to Sgr, as we discuss further in ii l7.1.5l 
Without any quantitative evidence for significant num- 
bers of Sgr stream stars that are close in both position 
and velocity to Segue 1, we do not see any reason to 
expect substantial Sgr contamination of our Segue 1 
member sample. 

7.1.4. Does the Sagittarius Stream Overlap m Velocity 
With Segue 1 ? 

NO09 also used SDSS observations of BHB stars to 
trace the kinematics of the observed Sgr debris near 
Segue 1 (while plenty of observations of Sgr stream ve- 
locities exist in other parts of the sky, the velocities near 
Segue 1 had not previously been measured) . They found 
that the main component of the stream has negative 
heliocentric (and galactocentric) velocities at this posi- 
tion, more than 200 km s~^ away from the velocity of 
Segue 1, as pointed out by G09. Another coherent BHB 
component, though, is present at much higher velocities, 
similar to the prediction from IFellhauer et al.l (2006) for 
older Sgr debris. NO09 concluded from this result that 
there is likely confusion between Segue 1 stars and Sgr 
stream stars in both position and velocity. However, 
even accepting for the moment that these stars are ac- 



2007t iKomivama et all 120071 : iSmolcic et al.ll2007[ ). but 

there is sti ff no c onsensus regarding their physical 
significance (jPenarrubia ct al..,2008„ ,2009.) . 



^^ While the ILaw fc Maiewskil I l2010al ) model provides a gener- 
ally reasonable match to Sgr stream data, the bifurcation of the 
stream in the SDSS foot print is not yet fu lly understood and is 
not present in the model. IPeiiarrubia et all (2010) argue that this 
structure can be a result of the original internal kinematics of the 
Sgr dwarf if its progenitor was a disk galaxy. 
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tually part of the Sgr stream (see ?? I7.1.3I and 17.1.51 for 
our counterarguments), two factors significantly dimin- 
ish this confusion. First, the BHB stream identified by 
NO09 does not actually appear to reach the location of 
Segue 1; it peaks at (5 ~ +24°, ~ 8° north of Segue 1, and 
seems to have petered out by the time it reaches Segue 1. 
Equally important is that the velocity of the BHB stars 
is Vhei ~ 195 km s~^ {Vgsr ~ 132 km s^^), offset from 
the velocity of Segue 1 by 13 km s^^ in the heliocentric 
frame and 19 km s^^ in the Galactic Standard of Rest 
(GSR) system. 

7.1.5. Contamination of the Segue 1 Member Sample 

Relying on the line of reasoning examined above, NO09 
proposed that Segue 1 is actually a star cluster whose 
derived properties have been distorted by contamination 
from Sgr stream stars. Such contamination is a difficult 
issue to quantify, because by definition any stars that 
could be contaminating the member sample must have 
very similar velocities, metallicities, distances, and ages 
to Segue 1 stars. The only way to assess definitively the 
expected number of contaminants would be with an even 
wider field survey to identify Segue 1-like stars that are 
far enough away from the galaxy so as to be very un- 
likely to be associated. While SDSS includes some of the 
desired data, the SDSS spectroscopic coverage of stars at 
faint magnitudes is very sparse, so the vast majority of 
stars do not have velocity measurements. Nevertheless, 
some do, and those observations can be used to estimate 
the significance of the contamination. 

The NO09 estimate of the contamination depends crit- 
ically on the assumptions discussed in ii l7.1.4l that the 
BHB stream they identify with Sagittarius is exactly 
coincident in both po sition and velocity with Segue 1. 
However, as noted in S I7.1.41 the surface density of BHB 
stars in the higher velocity component appears to be 
down by a factor of at least a few by the time it reaches 
Segue 1 (their Fig. 10), and the velocity offset compared 
to Segue 1 further reduces the contribution of these stars 
within the Segue 1 velocity selection window. 

To take into account these effects, we repeat the analy- 
sis described by NO09. Using DR7 data, if we select the 
BHB stars according to their heliocentric velocities, the 
positive velocity stream component has a mean velocity 
of Vfiei = 195 km s~^. Since the stream is extended spa- 
tially, and therefore likely has a velocity gradient along 
its length, it is more concentrated in the GSR velocity 
system, where its velocity is vgsr = 132 km s~^ and 
its velocity dispersion is 12.2 km s~^ (corresponding to 
an intrinsic dispersion of 10 km s~^ after the 7 km s~^ 
median velocity errors are removed). Given the velocity 
window spanned by the likely Segue 1 members (not in- 
cluding SDSSJ100704.35-I-160459.4) of 194.6 km s^^ < 
Vhei < 224.2 km s^^, only 41% of stars in the positive 
velocity stream would be expected to have velocities con- 
sistent with membership in Segue 1. Applying our ac- 
tual photometric selection criteria (§ 12. 2p instead of the 
broader selection box used by NO09 reduces our estimate 
of the surface density of Sgr stream stars at the declina- 
tion of Segue 1 to 120 deg~^ (for stars with r magnitudes 
between 17.5 and 21.7). The complete region of our spec- 
troscopic survey covers 0.087 deg^, and the effective area 
of the full survey (including the incompleteness at larger 



radii) is ~ 0.14 deg^. Assuming, as NO09 did, that half 
of these stars are in the negative velocity stream com- 
ponent, and removing the 59% of the positive velocity 
stream stars that would still lie outside the Segue 1 ve- 
locity range, suggests that a total of 2 — 3 Sgr stream 
stars could be in our sample. 

The next question is what effect including a few Sgr 
stars in an analysis of Segue 1 would have. We repeat 
the Monte Carlo simulation carried out by NO09 to an- 
swer this question. Using the same setup they did, with 
assumed velocity dispersions of 1 km s~^ for Segue 1 and 
10 km s~^ for the Sgr stream (and putting them both at 
the same mean velocity, contrary to the argument above) , 
we find that 5 Sgr stars must be included in the 71 star 
Segue 1 sample to have a significant chance (~ 20%) of 
boosting the apparent velocity dispersion of Segue 1 to 
at least 3.9 km s^^. Given the smaller number of con- 
taminating stars estimated above, we conclude that the 
inclusion of Sgr stream stars in the Segue 1 member sam- 
ple is not likely to provide the dominant component of 
the observed velocity dispersion. 

Also, when significant numbers of such contaminants 
are present they tend to have an easily visible effect on 
the velocity distribution (as they must if they are going 
to alter the dispersion). Visually, the simulated velocity 
histograms frequently appear to be composed of a narrow 
central peak containing most of the stars, surrounded by 
a few well-separated outliers (see Fig. |9]for an example). 
These outliers would raise suspicions in any membership 
classification scheme like the ones outlined in ii l3.1l and 
might well be discarded from the sample. Interlopers 
that happen to fall within the main peak of the velocity 
distribution (and are therefore more difficult to identify) 
do not have a significant impact on the velocity disper- 
sion; only stars in the wings of the distribution can both 
be mistaken for members and substantially change the 
apparent dispersion. However, our analysis shows that 
any such contaminating population cannot be large, and 
in Paper II we demonstrate that including an additional 
population does not change the derived parameters for 
Segue 1. 

Moreover, a closer examination of the positive ve- 
locity stream component calls into question the as- 
sumption that it is associated with Sagittarius at all. 
iNewberg et"all (|2010[ ) used BHB stars in SDSS and the 
SEGUE survey to trace the Orphan Stream across the 
sky and noted that it passes slightly north of Segue 1, 
at a distance of ~ 25 kpc and a velocity of vgsr = 
130 km s"^ Isolating the SDSS BHB stars within 1 a of 
the positive velocity stream's mean velocity, we find that 
their spatial distribution closely matches the track of the 
Orp han Stream det e rmine d by [Belokurov ct al. (2007q) 
and INewberg et al.l (|2010( ) , as shown in Figure [ini The 
good correspondence between the path of the Orphan 
Stream and the BHB stars at the same velocity is highly 
suggestive that these stars are members of the Orphan 
Stream rather than an old (and heretofore undetected) 
wrap of the Sagittarius stream. If this peak is indeed 
related to the Orphan Stream, which is narrow and spa- 
tially confined, and not Sagittarius, then it is quite un- 
likely that any main sequence stars associated with this 
feature would be located close enough to Segue 1 to ap- 
pear as contaminants in our survey. That would then 
imply that most or all of the Sgr stars near Segue 1 are at 
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Fig. 9. — Monte Carlo simulation of the expected velocity distri- 
bution that would be obtained in the presence of significant con- 
tamination by the Sgr stream. We assume intrinsic dispersions of 
1 km s~^ for Segue 1, 10 km s^^ for the Sgr stream, mean veloci- 
ties of 208 km s~^ for both components, median velocity errors of 
5 km s~^ (with a minimum of 2.2 km s~^), and 5 Sgr stars in a 71 
star sample. The velocities of the full simulated sample are shown 
in black, with the Sgr contaminants overplotted in cyan and the 
observed Segue 1 stars in red. 



negative velocities, as suggested in §[7T|3l which would 
further reduce our estimate of the possible contamination 
by Sgr above. 

After considering each piece of evidence in concert, we 
thus conclude that contamination by the Sgr stream does 
not provide a very plausible explanation for the large 
velocity dispersion of Segue 1. 

7.2. Signatures of Tidal Disruption 

Having determined the nature of Segue 1 , the effect of 
binary stars on the velocity dispersion, and the level of 
contamination by the Sgr stream, the final issue we must 
analyze is the impact of Milky Way tides. Unfortunately, 
while the presence of tidal tails would be incontrovertible 
evidence of tidal effects, the contrapositive is not true: 
there are no observations that can conclusively rule out 
tidal disruption. We therefore consider several possible 
signatures of tides. 

• First, our spectroscopic survey shows there are 
no obvious tidal tails connected to Segue 1 (see 
Figure [3]). Although the spatial distribution is 
not uniform, we find Segue 1 members in ev- 
ery direction around the galaxy rather than the 
bipolar pattern that tidal tails would be ex- 
pected to produce. The apparent dumpiness of 
the member stars toward the western edge of 
the galaxy may si mply be the result of small 
number statistics (jMartin et al.l HOOS). While 
the spectroscopic member sample confirms that 
Segue 1 has an elliptical shape, nonzero elliptic- 
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Fig. 10. — Spatial distribution of the BHB stars that can be 
confidently associated with the positive velocity stream identified 
by NO09. The green star marks the position of Segue 1, and the 
cyan and blue l i nes in dicat e the traces of the Orphan Stream from 
INewberg et al.l l|20ia ) and IBelokurov et al.l (|2007bl 'l , respectively. 
The close correspondence between the path of the Orphan Stream 
and the positions of the BHB stars at the same velocity suggests 
that these stars are members of the Orphan Stream rather than 
the Sagittarius stream. 

ities arc not necessarily associated with tidal in- 
fluences (Mmioz, Maicwski, & Johnston 2008) and 
may just reflect the shape with which the galaxy 
formed. We also note that the tidal tails seen in the 
SPSS photometry in th e Segue 1 discovery paper 
(JBelokurov et al.ll2007al) have not been co nfirmed 
by deeper followup (JBelokurov et aIll2007aL NO09, 
Mufioz et al., in prep.). 

• Velocity gradients are a commonly used indicator 
of tidal disruption, although like tidal tails they 
may only be visible in particular g e ometries and 
at lar ge radii (iPiatek fc Prvoii 119951 : iMunoz et al.l 
120081 : iLokas et al.ll2008[ ). We see no velocity gradi- 
ent across Segue 1 ; the mean velocities of the stars 
in the eastern and western halves of the galaxy 
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agree within their uncertainties. The inclusion or 
exclusion of the ambiguous members, RR Lyraes, 
and binaries does not affect this result. There 
are also no apparent trends in the mean velocity 
from one end of the gala xy to the o ther (unlike, 
e.g., WiUman 1; Willman et al.ll2010D. Using the 
methodology described bv lStrigaril (|201Cil ) to calcu- 
late completely general constraints on the rotation 
of Segue 1, we derive a 90% confidence upper limit 
of 5.0 km s~^ on the rotation amplitude, but we 
note that samples a factor of 2 larger than ours are 
typically required to detect rotation. 

• A velocity dispersion profile that rises at large radii 
is often regarded as a possible result of tidal strip- 
ping, although the same behavior can be inter- 
preted as evidence for an extended dark matter 
halo as well. Our member sample is not large 
enough to divide the data into more than a few 
radial bins, but we can begin to investigate the 
shape of the dispersion profile. Velocity dispersion 
profiles with two different binnings are displayed 
in Figure [TTJ As seen in the left panel of Fig- 
ure m the velocity dispersion appears to reach a 
local minimum at a radius of ~ 3' before increas- 
ing back to its central value at larger radii. This 
shape seems to be independent of the exact binning 
chosen (and in fact is visible in the unbinned data), 
but as the error bars in the figure show, it is not 
statistically significant. While it is possible that 
the dispersion increase could suggest that the stars 
beyond ~ 8' from the center of the galaxy have 
been stripped, we caution that apparent features in 
other data sets of similar size (or even larger) have 
often disappeared when larger sa mples of velocity 
measurements become a vailable (| Wilkinson et al.l 
120041: iKlevna et aIll2004D . With the modest num- 
ber of bins possible for a sample of 71 members, we 
view the shape of the dispersion profile of Segue 1 
as possibly interesting but not necessarily mean- 
ingful at this point. It is also worth pointing out 
that the mass implied by the central velocity dis- 
persion, even if the decline at ^ 3' is real, is enough 
to put the tidal radius beyond the observed extent 
of the galaxy (see below), suggesting a consistency 
problem for the tidal interpretation. 

• Finally, "extratidal" excesses of stars at large radii 
are frequently con sidered to b e indicative of tidal 
disturbances. The iKin^ (|1962t ) tidal radius or lim- 
iting radius of Segue 1 is not well known because of 
the lack of deep enough wide-field photometry (al- 
though NO09 estimate a value of ^ 26'), but our 
complete spectroscopic sample enables us to inves- 
tigate the stellar profile out to r ~ 13'. Our obser- 
vations are effectively complete within 2 half-light 
radii of the center of the galaxy fJi l2.6l) . where we 
identify 61 member stars. Between 2 and 3 times 
the half-light radius we obtained successful spectra 
for 62 out of 75 stars located in the highest pri- 
ority photometric selection region, for a complete- 
ness of 83%. We therefore adjust the 9 observed 
members in that annulus to a projected total of 11 
members, yielding 72 member stars within 3 half- 
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Fig. 11. — Velocity dispersion profile of Segue 1. The open 
squares show the profile obtained for bins of ~ 23 stars each, and 
the filled red circles show the profile for bins of ~ 15 stars. While 
the decrease in the velocity dispersion at intermediate radii (also 
visible in Figure^ does not appear to be an artifact of the binning, 
it is only significant at the 1 a level. 

light radii (13.2'). We find 39 member stars within 
1 ''half (54% of the total), 22 member stars between 
1 ''half and 2 rhaif (31%), and estimate 11 member 
stars between 2 rhaif and 3 rhaif (15%), compared 
to the expected numbers of 56%, 33%, and 11% for 
a Plummer profile (once the 10% of the stars that 
should lie beyond 3 rhaif are removed from consid- 
eration). The radial profile of Segue 1 thus does 
not show any excess out to at least 3 rhaif (88 pc). 

The above arguments demonstrate an absence of evi- 
dence in favor of tidal disruption, but as mentioned at 
the beginning of th is section and discussed in detail by 
IMuhoz et al.l (J2008D . none of them (singly or in concert) 
are sufficient to prove that Segue 1 is not being tidally 
disrupted. Perhaps the strongest evidence for the ab- 
sence of tidal effects results from consistency checks be- 
tween the mass we measure for Segue 1, the correspond- 
ing tidal radius, and the timescale for tidal disruption. 

As discussed by G09, the current position of Segue 1 
is difficult to reconcile with a scenario in which it is in 
the final stages of disruption. Segue 1 has a crossing 
time of ~ 10 yr, and at a velocity of ~ 200 km s~^ it 
will travel less than 2 kpc per crossing time. In order 
for the observed kinematics of Segue 1 to be significantly 
distorted by tides, the galaxy must be within a few cross- 
ing times of its orbital pericenter. Conservatively, then. 
Segue 1 should be no more than ~ 10 kpc past pericen- 
ter. Segue 1 is located 28 kpc from the Galactic Center, 
which would place its pericenter at a Galactocentric dis- 
tance of at least 18 kpc, inconsistent with the closest 
approach to the Milky Way that would be required to 
disrupt it (see below). 

For the lAU value of the Milky Way rotation veloc- 
ity (220 km s~^), the mass enclosed at the position of 
Segue 1 is 3 X 10^^ M0. The mass of Segue 1 is b est con- 
strained at the half-light radius of the galaxy ( Wolf et al.l 
[2OT0l) . where we obtain Afhait = ^-^tH x 10^ Mq. 
Even if we assume (without any physical basis) that 
the mass distribution is arbitrarily truncated at the 
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half-light radius, the inst antaneous Jacobi radius (e.g., 
iBinnev fc Tremaind [20081 Equation 8.91) for Segue 1 
would be ~ 250 pc. If we use the Jacobi radius as an esti- 
mate of the tidal radius, then all of the stars we observed 
are well within the present-day tidal radius. Although 
IBinnev fc Tremaind present a detailed discussion of why 
the Jacobi radius is necessarily an imperfect estimator, it 
is worth noting that subhalos in the Aquarius simulation 
(jSpringel et al.|[2008() show a strong correlation between 
the radius of the subhalo and the Jacobi radius, sup- 
porting the use of the Jacobi radius as the tidal radius in 
practice. In order to bring the tidal radius in to the posi- 
tion of our outermost confirmed member, the pericenter 
of Segue I's orbit must be less than ~ 10 kpc. Substan- 
tially disturbing stars at the half-light radius, where the 
mass is being measured, requires an orbital pericenter of 
less than ~ 4 kpc (eccentricity greater than 0.75, since 
Segue 1 is clearly not near its apocenter at the present 

time). The Jacobi radius scales as Mg'^ ^, so even ma- 
jor revisions to the derived mass do not affect our con- 
clusion. This calculation is quite conservative, because 
it relies most strongly on the central kinematics of the 
galaxy where the observational uncertainties are smallest 
and tidal effects are weakest. 

Incorporating reasonable assumptions about the extent 
of the dark matter halo of Segue 1 only strengthens this 
result. In the Via Lactea II simulation ([Diemand et al.l 
|2008[) . subhalos with Knax > 10 km s'^ (see §[8]) that cur- 
rently reside between 20 and 40 kpc from the host halo 
have median tidal truncation radii of ^ 500 pc. Since 
these simulations self-consistently include tides and or- 
bital trajectories, there is good reason to suspect that 
the mass of Segue 1 extends well past rhaif- If we ex- 
trapolate the Segue 1 mass beyond the observed region 
using CDM priors (since current simulations cannot re- 
solve radii smaller than ^ 100 pc), we find a mass within 
100 pc of Afioo = 2.2 X 10^ Mq and a mass within 300 pc 
of M3Q0 = 1.4 X 10^ Mq, consist ent with the commo n 
mass scale of Milky Way satellites (jStrigari et al.ll20d8al ). 
With these larger masses, the current tidal radius would 
increase to 400 — 700 pc, making the center of Segue 1 
nearly impervious to tides for any plausible orbit. 

8. THE IMPORTANCE OF SEGUE 1 FOR DARK 
MATTER STUDIES 

The large estimated mass of Segue 1 and its very 
small size (it has the smallest half-light radius of any 
known Local Group dwarf galaxy, with the possible 
exception of Willman 1) mean that Segue 1 also has 
the densest known concentration of dark matter. The 
average density enclosed within its half-light radius is 
2.5l^g M© pc""^, substantially higher than that found in 
other dwarf galaxie s (IGilm ore et al.ll2007tlSimqn fc Gehal 
[2OO7I: IWalker et al.1 !2009a: ToUerud et apl2011l [Fig. 17n 
or the solar neighborhood fe.g.. lBahcalllll98^ r For com- 
parison purposes, this density is equal to the ambient 
density of dark matter at z ~ 300 and the average density 
of objects that collapsed at z ~ 50. Given the extremely 
high mass-to-light ratio of Segue 1 (§ 15. Sp . it is safe to 
equate the dark matter density with the total density. 

In the context of ACDM, densities as high as those 
that we infer within the half-light radius of Segue 1 
are indicative of massive subhalos. Using the mass es- 



timator from iWolf et ahl ()2010D . the circular velocity 
at the half-light radius is related to the measured ve- 
locity dispersion via F(rhaif) = VScr ~ 6.4 km s~^. 
This demands that T^nax > 6.4 km s~^ for the halo 
hosting Segue 1. Convolving the central circular ve- 
locity with CD M-based priors sug gests that 14iax > 
10 km s-i (;e.g.. lBullocket al.ll2010l ). Such halos are in- 
deed found in the mo st advanced current N-body simula- 
tions (Die mand et a l. 2008; Springel et al. 2008), demon- 
strating that the density of Segue 1 is reasonable in a 
ACDM universe, but how common it is for galaxies of 
Segue I's luminosity to be found in very massive subha- 
los is not yet clear. 

Since the flux of high-energy particles from dark mat- 
ter annihilation scales as p^-^r^ /d? and Segue 1 is also 
the second-nearest dwarf galaxy to the Sun, Segue 1 
is clearly a hig h priority target for indirect detec- 
tion experiments (l^fartinez et ai.ll2009l: lEssig et al.|[200l 
iScott et al.ll2010f l lEssig et al.l (|20l'or use the sample of 
member stars presented here to carry out more detailed 
calculations of the expected gamma-ray and neutrino flux 
from dark matter annihilation in Segue 1. That analy- 
sis shows that Segue 1 is expected to be among the two 
brightest sources of annihilation radiation from Milky 
Way satellites, and may be the brightest known dwarf 
galaxy. We strongly encourage future indirect detection 
searches for dark matter to target Segue 1. 

Finally, the high density of Segue 1 provides impor- 
tant leverage for constraints on the phase-space density 
of dark matter particles, which is often estimated by the 
related quantity QdmIt) = puuir) / (J-Duir)^ defined by 
iHogan fc DalcantonI ()2000l ). Unfortunately, Qdm can- 
not be measured directly from velocity dispersion data. 
While pdm may be determined fairly accurately within 
'"half, cp M is not observab le. Generally, we expect ctdm > 
a^, (e.g. JWolf et al.|[2010[ ) because the dark matter veloc- 
ity dispersion is governed by the total mass beyond the 
stellar radius (which cannot be measured). This implies 
that pdm/c* provides only an upper limit on Qdm at 
any particular radius. Some caution is advisable when 
reading the literature on this subject. 

Given that Qdm cannot be measured directly, we must 
rely on model fitting inspired by a theory prior in order 
to quantify Qdm constraints from Segue 1. Of particular 
interest is the case of warm dark matter (WDM), where 
the primordial phase-space density may produce obser- 
vationally accessible cores in the dark matter density. 
It is therefore useful to assume a dark matter density 
profile that is compati ble with WDM rather than the 
usual CDM Navarro, F renk. fc Whitd (fT99g) or similar 
profile. We use a cored isothermal profile for illustra- 
tion, which should have approximately the right shape 
for such models. With this profile and the kinematic 
data presented in this paper, we determine the posterior 
probability density for Qdm at equally spaced logarith- 
mic radii out to the stellar tidal radius following, e.g., 
iStrigari et a l.l (I2008B)- We assume uniform priors on 
the scale radius and scale density for the isothermal pro- 
file, and we additionally make the assumption of velocity 
isotropy for the stars and the dark matter. Under these 
assumptions, we find a lower limit on the central value 
for Qbm to be ~ 10~^ M© pc~^ (km s~^)~^, higher than 



Complete Spectroscopic Survey of Segue 1 



19 



that determined for any other galaxvF^ Estimates of this 
kind can place a lower limit on the allowed mass range 
for various WDM candidates, and we suggest that more 
detailed treatments of the phase-space density in Segue 1 
to quantify these constraints would be very worthwhile. 

9. SUMMARY AND CONCLUSIONS 

We have presented a comprehensive spectroscopic sur- 
vey of the ultra- faint Milky Way satellite Segue 1. The 
observations were designed both to search for potential 
tidal debris around Segue 1 and to constrain the effects 
of binary stars on its velocity dispersion. Within a radius 
of 10' (67 pc) from the center of Segue 1, we measured 
the velocities of 98.2% of the candidate member stars. 
We identified 71 likely members, which we used to study 
the metallicity, kinematics, and nature of Segue 1. 

The six red giants in Segue 1 have a mean metallicity 
of [Fe/H] = —2.5, and span a range of nearly 2 dex from 
[Fe/H] = -3.4 to [Fe/H] = -1.6. Both the presence 
of extremely metal-poor stars with [Fe/H] < — 3 and the 
enormous metallicity spread demonstrate unambiguously 
that Segue 1 is a galaxy, rather than a globular cluster 
as some previous studies have suggested. 

Using the ne w Bayesian method pr esented in our com- 
panion paper (jMartinez et al.l 120101 ) , we analyzed the 
kinematics of the entire observed data set, allowing for 
contamination by Milky Way foreground stars and em- 
ploying repeated velocity measurements for a subsample 
of the targets to correct for the effect of binary stars. We 
derived an intrinsic velocity dispersion of 3.7_2]^ km s~^ 
for Segue 1 and only a 2% probability that the dis- 
persion is small enough to be provided by the stellar 
mass of Segue 1 alone. The estimated mass contained 
within the half-light radius is 5.8^31 x 10^ Mq, giving 
Segue 1 a V-band mass-to-light ratio at that radius of 
~ 3400 Mq/Lq. 

Based on updated data and models, we re-examined 
earlier proposals that Segue 1 is tidally disrupting and 
that kinematic studies of it are likely to be contaminated 
by the Sagittarius stream. We showed that there is no 
observational evidence supporting the possibility of tidal 
disruption, and that the tidal radius of Segue 1 has likely 
always exceeded its stellar extent unless it has an orbital 
pericenter around the Milky Way of less than ~ 4 kpc. 
We also determined that contamination by Sgr stream 



stars is significantly lower than previously estimated; our 
current member sample is unlikely to contain more than 
3 contaminants, which is not enough to substantially in- 
fiate the velocity dispersion. 

Taken together, the results of our observations clearly 
point to the interpretation that Segue 1 is a dark matter- 
dominated galaxy — in fact, it has the highest mass-to- 
light ratio, and is therefore the darkest galaxy, yet found. 
The mean density inferred for Segue 1 within its half-light 
radius is consistent with the extrapolated density pro- 
files of massive subha los in high resolution ACDM galac- 
tic halo simulations (jMadau. Diemand. fc KuhJenl 120081 : 
iSpringel et al.|[2008l) . The relative proximity of Segue 1 
makes a strong case for considering Segue 1 in future 
searches for the products of da rk matter annihilation pro- 
cesses (e.g.. lEssig et al.ir2010r ). The density of dark mat- 
ter within the inner 38 pc of Segue 1, 2.5l^'g M© pc^"^ or 
~ 100 GeV/c^ cm~'^, is the highest dark matter density 
yet determined, and consequently has broad implications 
for particle physics models and galaxy formation on small 
scales. 
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TABLE 3 

Segue 1 Velocity Measurements 



Star 


Velocity 


SCa 


Radius 


MJD 


g 


r 


i 


Member 


EM Member 


Bayesian 




(kms-i) 


(A) 


(arcmin) 










(Subjective)^ 


Prob.b 


Member Pre 


SDSSJ100613. 98+155436.1 


-46.8 ±4.6 


3.3 ±0.6 


15.4 


54890.243 


24.28 


22.64 


21.17 





-9.999 


-9.999 


SDSSJ100614.24+160424.7 


-66.7 ±2.4 


2.4 ±0.3 


11.8 


54890.243 


20.29 


19.97 


19.81 





0.000 


0.000 


SDSSJ100614.31+160050.9 


150.1 ±3.3 


2.9 ±0.5 


12.3 


54890.243 


21.52 


21.03 


21.13 





0.000 


0.000 


SDSSJ100614.40+160013.0 


70.1 ±2.2 


4.8 ±0.3 


12.5 


54890.243 


17.14 


16.44 


16.16 





0.000 


0.000 


SDSSJ100614.78+155512.7 


45.2 ±2.9 


1.5 ±0.4 


14.8 


54890.243 


24.83 


22.50 


21.06 





-9.999 


-9.999 


SDSSJ100614.87+160858.7 


29.1 ±2.9 


4.4 ±0.4 


12.5 


54890.243 


20.09 


19.55 


19.32 





0.000 


0.000 


SDSSJ100615.16+155556.6 


54.3 ±2.6 


2.9 ±0.3 


14.3 


54890.243 


22.30 


21.07 


19.96 





-9.999 


-9.999 


SDSSJ100615.53+160056.9 


-6.7 ±2.7 


1.7 ±0.4 


12.0 


54890.243 


21.04 


20.73 


20.65 





0.000 


0.000 


SDSSJ100616.95+160524.3 


-10.0 ±2.4 


3.3 ±0.3 


11.2 


54890.243 


23.22 


21.47 


20.39 





-9.999 


-9.999 


SDSSJ100617.35+155606.6 


1.7 ±2.4 


2.1 ±0.4 


13.8 


54890.243 


22.81 


21.47 


20.68 





-9.999 


-9.999 



Note. — This table will be published in its entirety in the electronic edition of the Journal. A portion is shown here for guidance regarding form ; 
content. 

'^ Member status according to the criteria established at the beginning of § 13.11 1 indicates membership, and is for non-members. 

'' Membership probability from the EM algorithm. The algorithm is run on the subset of stars whose colors and magnitudes are consistent with mcmbcrsl 
so photometric non-members are indicated by probability —9.999. 

'^ Membership probability from the Bayesian approach. As with the EM algorithm, these calculations are run on the subset of stars whose colors ; 
magnitudes arc consistent with membership, so photometric non-members are indicated by probability —9.999. 



Complete Speetroscopic Survey of Segue 1 



23 



TABLE 4 

Segue 1 Metallicity Measurements 



Star 



Tnf 



logs 



[Fc/H] 



Number of 
measurements 



Dispersion between 
measurements 



SDSSJ100702.46+155055.3 
SDSSJ100714.58+160154.5 
SDSSJ100652.33+160235.8 
SDSSJ100742.72+160106.9 
SDSSJ100710.08+160623.9 
SDSSJ100639.33+160008.9 



18.48 


17.94 


17.73 


5148 ± 102 


2.64 


-2.48 ±0.15 


2 


0.42 dex 


18.83 


18.30 


18.06 


5102 ± 109 


2.78 


-1.73 ±0.14 


3 


0.07 dex 


18.87 


18.39 


18.16 


5271 ± 132 


2.85 


-3.40 ±0.17 


3 


0.32 dex 


19.72 


18.59 


18.14 


5251 ±111 


2.76 


-2.50 ±0.14 


1 




19.20 


18.71 


18.45 


5106 ± 109 


2.96 


-1.63 ±0.14 


4 


0.03 dex 


19.44 


19.03 


18.90 


5643 ± 186 


3.21 


< -3.4 


3 
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APPENDIX 

APPENDIX: COMMENTS ON THE MEMBERSHIP OR NON-MEMBERSHIP OF INDIVIDUAL 

STARS 

In this appendix, we briefly discuss the case for membership or non-membership for several stars whose membership 
status is not clear-cut. Unlike SDSSJ100704. 35-1-160459. 4, however, the decision of whether to include or exclude 
these stars has no appreciable effect on the derived properties of Segue 1. We remind the reader that in the Bayesian 
approach used for our main results all of these stars (except the last two) are included in the calculations, weighted 
according to their membership probabilities. Other schemes generally require each star to be classified as either a 
member or non-member. 

• SDSSJ100637.49+161155.1 — Very much like SDSSJ100704.35-hl60459.4, SDSSJ100637.49-I-161155.1 sits per- 
fectly on the Segue 1 fiducial sequence but has a significantly higher velocity (more than 3cr even for the largest 
possible value for the dispersion) than the galaxy. Given its relatively large distance from the center of Segue 1 
(9.7') and its higher than average CaT EW {W = 3.66 A), we consider this star unlikely to b e a member. Its 
EM membership probability is 0.69 and its membership probability from [Martinez et al.l ()2010l ) is 0.23. 

• SDSSJ100622.85+155643.0 — SDSSJ100622.85-hl55643.0 has a very large uncertainty on its velocity measure- 
ment {v = 223.9 ± 37.8 km s~^) because of the low S/N of its spectrum and an apparently broad Ha line. Since 
much of the Icr range for its velocity puts the star at a velocity consistent with Segue 1, we subjectively classi- 
fied it as a member, but the statistical algorithms give it lower probabilities (p = 0.73 for EM and p = 0.50 for 
[Martinez et al.|[2010l However, the large velocity error means that SDSSJ100622. 85-1-155643.0 receives essentially 
no weight in determining the Segue 1 velocity dispersion, so its true membership status is not important. 

• SDSSJ100711. 80+1606304 — SDSSJ100711.80-H160630.4 lies well outside the Segue 1 velocity range at u = 
247.1 km s^^, but also has a large velocity uncertainty of 15.9 km s~^. The very high velocity led us to classify it 
as a non-member, but the objective techniques recognize that there is a reasonable chance that the star's actual 
velocity could be significantly lower, which might make it a member. Weighting its position near t he center of 
Segue 1 quite heavily, the EM algorithm gives a membership probability oi p — 0.98, while the Ma rtinez et al.l 
algorithm more conservatively estimates p = 0.70. Again, the large velocity error minimizes its impact on the 
derived velocity dispersion. 

• SDSSJIOO743.55+I6O947.2 — SDSSJ100743. 55-^160947.2 has a velocity near the high end of the Segue 1 velocity 
range {v = 223.4 ± 5.3 km s"""^), but the uncertainty is large enough for it plausibly to be a member. In the 
color-magnitude diagram, it is located just outside the highest priority selection region, but again close enough 
that the photometric errors certainly allow it to be a member. We consider this star to be a probable member, 
but the objective algorithms give it moderate membership probabilities (p = 0.87 for EM and p = 0.54 for the 
Bayesian determination) because of its high velocity and large radius. 

• SDSSJ100630.96+155543.6 — SDSSJ100630.96-M55543.6 has a velocity several standard deviations smaller 
than the systemic velocity of Segue 1, although the large velocity uncertainties (8.7 and 11.8 km s~^ for the two 
measurements) make this differen ce of marginal significance. Combined with its position 11.6' from the center 
of the galaxy, the [Martinez et all algorithm hedges its bets at p = 0.60, while the EM membership probability is 
0.94. This star is another probable member. 

• SDSSJ100658. 11+160701. 4 — With a velocity of 187.0 ± 2.3 km s^^ and a 5 - i color - 0.25 mag blueward of 
the red giant branch (outside both the high and low priority CMD selection regions), SDSSJ100658. 11-1-160701. 4 
appears to be a clear non-member star. However, the star could conceivably be an SX Phoenicis variable, in 
which case membership in Segue 1 would be possible, although as a variable star it would still be excluded from 
our analysis of the kin ematics. SD SSJ1006 58. 11-1-16 0701.4 is redder than would generally be expected for an 
SX Phoenicis star fe.g.. l01ech et al...2005. : .Moretti et al..,2009.) . but we cannot rule out such a classification with 
the available data. 



• 



SDSSJ 100700. 75+ 160300. 5 — In our standard reduction of the data, using the arc spectrum obtained closest to 
the time of the observations, we were not able to determine a reliable velocity for SDSSJ100700. 75-1-160300. 5. The 
cross-correlation with the best-fitting template spectrum produced multiple widely separated cross-correlation 
peaks, with no obvious way to identify the correct solution. Switching to the arc frame from the end of the 
night (^ 10 hours after the mask was observed) instead of the one from the afternoon [^ 3 hours before the 
observations) produced a cleaner spectrum with a velocity of u = 208.4 ± 9.7 km s""'^. Because none of the 
other spectra on this mask required such special treatment, and the increased time between calibrations and 
observations provides more opportunity for changes in the instrument, we regard this measurement as somewhat 
questionable and omit the star from our sample. It is probably a member of Segue 1, but with a velocity 
consistent with the systemic velocity of the galaxy and a large velocity error, including it would not change any 
of our results. 



